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Introduction 

A FORTRAN computer program has been developed which automatically 
targets two and three burn rendezvous missions and performs feedback guidance 
using the previously developed GUIDE algorithm. The program was designed to 
accept a large class of orbit specifications and automatically chooses a two 
or three burn mission depending upon the time alignment of the vehicle and 
target. The orbits may be specified as any combination of circular and ellip- 
tical orbits and may be coplanar or inclined, but must be aligned coaxially 
(i.e. line of intersection of orbital planes and orbital major axes coincident) 
with their perigees in the same direction. The program accomplishes the re- 
quired targetting by repeatedly converging successively more complex missions. 

It solves the coplanar impulsive version of the mission, then the finite burn 
coplanar mission and finally the full plane change mission. The GUIDE algorithm 
is exercised in a feedback guidance mode by taking the targeted solution and 
moving the vehicle state step by step ahead in time adding acceleration and 
navigational errors and reconverging from the perturbed states at fixed guid- 
ance update intervals. 

The targetting and guidance algorithm converges all two burn missions 
easily and exhibits good guidance behavior for these missions. Three burn mis- 
sions were much more sensitive and required special loops to insure convergence'. 
The outbound three burn mission had to be converged backwards in time and plane 
change was most readily incorporated by eliminating the third burn and solving 
the appropriate two burn mission, reintroducing the third burn at the end. In 
a targetting mode these techniques cause no particular problem and insure con- 
vergence. In guidance mode the convergence problems are more difficult to com- 
pensate for and may limit real time use. The program as it now stands attempts 
to optimize over all three burns and although it has maintained convergence for 
all missions attempted, the guidance corrections have been larger than desired. 
In the future it may be necessary to solve the guidance problem over the first 
burn as a rendezvous with the desired phasing or transfer orbit and to only 
introduce the third burn after completion of the first one. 

Another study that needs to be undertaken is to optimize the soft con- 
straint weights using the Monte Carlo capability built into the. program. By 
altering the weights and noting the tradeoffs made between burn time and orbital 
injection error, a better estimate of optimal soft constraint weights can be 
obtained. - 
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The remainder of this document describes the targetting and guidance 
program in detail, giving an overview of the program control and organization, 
a summary of program inputs and outputs and a detailed description of each of 
the subparts of the program. Also included in the document is a description 
of the GUIDE subroutine BVAL5, which was, .altered to incorporate the sofi con- 
straint formulation, and is fully documented. The other GUIDE subroutines are 
essentially the same as the ones described in the GUIDE 71/6 document^ and 
are not described here. 



^ Cohen, A.O., "Guide 71/6 Program Documentation", IBM Federal Systems 
Division, Burlington, Mass., October 4, 1971. 
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Program Overview 

The program is controlled by routine MAIN, which oversees the impulsive 
targetting, the convergence of the orbital transfer, and the feedback guidance. 
The impulsive targetting is accomplished by first determining the elements of 
both orbits, then defining the transfer orbit and phasing orbit (3 burn only) 
and dete»mn‘ning the velocities at apogee and perigee of each orbit, rlexc the 
delta v's are calculated and the burn and coast times calculated. The transfer 
orbit is chosen to be tangent at both end points to the principal orbits, and 
Jthe_missicn is classified as inbound or outbound depending on whether apogee 
of the final (target) orbit is less than cr greater than apogee of -the initial 
jorbit. The phasing orbit is chosen to lie as close as possible to the one which 
results from splitting the burn at perigee into two equal halves. A closed- 
form solution is used for initial costate. 

The converged finite-burn solution is arrived at by repeatedly conver- 
ging successively more complex missions, starting with a planar mission and 
gradually adding in the plane change required (10° steps). To maintain con- 
vergence for outbound 3-burn missions, it was necessary to rearrange each mis- 
sion and converge it in a backwards fashion, from the target orbit to the 
vehicle (initial) orbit. The plane change mentioned above was facilitated by 
changing the 3-burn mission to a 2-burn mission where the planar-converged 
phasing orbit was substituted for the closer orbit. After converging the 2- 
burn mission with the total plane change, the 3-burn mission was reinstated 
and converged. Finally, the 3-burn outbound mission is turned around to its 
normal mode and reconverged. 

After targetting has been done, the guidance portion of the program 
is run in a feedback mode, in which it is made to respond to simulated pertur- 
vations. The routine MAIN calls BCBCB or CSCB to propagate the vehicle along 
each arc of the mission, and Monte Carlo statistics are collected at appropri- 
ate points and summarized at the end. 

Further details of the operation of the program, as well as the routines 
employed, are described in the pages which follow. 
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User's Guide 

The program is set up using NAMELIST input for ease of operation. 

This allows default parameter values to be specified and reduces the amount 
of input necessary for program execution. Typical space tug vehicle para- 
meters are hard coded as default values and tug missions can be performed by 
simply specifying the desired initial and final orbits. The basic program 

,i'i 

philosophy is to use the. orbital definitions to define whether the mission 
will be two or three burns. If the mission is circular to circular coplanar* 
or if the orbital elements are defined with no positions along the orbits given, 
or if the positions of the vehicle and target allow a two burn rendezvous, a 
two burn orbital transfer will be defined. Under all other conditions three 
burn transfers will be used. The integer NOTARG is used to control which por- 
tions of the program are executed. If N0TARG=-1 only targetting is performed. 

If N0TARG=1 a converged solution for the orbital transfer is read in using 
NAMELIST NAMSL2 and only the feedback guidance part of the program will be 
executed. If NOTARG is any other value both the targetting and guidance will 
be performed. The inputs and outputs and individual subroutines will be des- 
cribed in detail in the sections which follow. 

Program Inputs 

The- program ■inputs~are’broken”tnto“‘three"basic“groups':'”those"whi ch'^de''^ 
fine the vehicle's capabilities, those used to specify the initial and final 
orbits, and those used to define the Monte Carlo and perturbation parameters 
needed for feedback guidance evaluation. 

A. Vehicle Constants 

The following parameters are used to specify the vehicle, and must be 
in metric units. If specific impulse is inputted it Is used to calculate mass 
rate. The default values for the parameters are typical of a space tug con- 


figuration. 




Name 

Symbol 

Definition 

Default Value 

AM0 

^0 

Initial vehicle mass in kg 

28803.1155 kg 
.. (63500 lbs) 

THRUST 

T 

Thrust in kilo-Newtons 

66.7233 kn 
(15000 lbs) * 
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Name Symbol Definition . Default Value 

Specific impulse in seconds 440 sec 
m Mass rate in kg/sec 15.4634 kg/sec 

Orbit Specifications 

The vehicle and target orbits may be specified in four separate v^fays 
listed as sets 1-4 below. (It is assumed that both will be specified in the 
same fashion.) For all of the orbital definitions the perigee directions must 
be equal and coincident with the line of intersection of the orbital planes. 

If sets 2, 3 or 4 are used to specify the orbits, these conditions are satis- 
fied automatically due to the way the orbital positions and coordinate systems 
are defined. If position and velocity vectors and times (set 1) are specified, 
the program will test to see that the conditions are satisfied and will stop 
if the proper perigee and line of nodes alignment is not found. When set 1 is 
used to specify the data the relative inclination between orbits is measured 
from vehicle to target orbit at perigee. In all other cases relative inclin- 
ation is set by the input data. If sets 2, 3 or 4 are used to specify the or- 
bits and the true anomalies (TANOM0 and TANOMT) are greater than or equal to zero, 
they will be used to specify the orbital positions. If true anomalies are not 
-specified and T0- and' TT are"greater"than"or“eqtraT“to—zero“ttrey“wrH“be-assumed 
to be mean anomalies and used to specify the orbital positions. If neither of 
the anomalies are specified, the orbital positions will be arbitrarily chosen 
to allow a two burn rendezvous. If no complete set of vehicle and target orbital 
data is available, the program will print the existing data and stop. 


Name 

Symbol 

Definition 

Default Value 

”r0(3),V0{3)J0 


Position (km) and velocity 
(km/sec) vectors at t^ for 
vehicle orbit 

T0=-1 

_ RT(3), VT(3), TT 


Position (km) and velocity 
(km/sec) vectors at t, for 
target orbit 

— 1 
li 


SPFIMP 

AMDOT 
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Name 

Symbol 

Definition 

6 

Default Value 


A0, E0 

aQs Sq 

Semi-major axis (km) and eccen- 
tricity for vehicle orbit 

0, -1 

■ 

s 

f- 

AT, ET 

- 


Semi-major axis (km) and eccen- 
tricity for target orbit 

0, -1 

E 

T 

2 

RELINC 

i 

Signed relative inclination (deg) 
as measured from 'vehicle to target 
orbit 

0 


TANOM0JANOMT* 

f 

True anomalies (not required) (deg) 

-1 


T0, TT * 

M 

Mean anomalies if true anomalies 
not specified (not required) (sec) 

-1 


* described more fully 

in text above 



HAP0, HPG0 

' 

- 

Height at apogee and perigee for 
vehicle orbit (km) 

None 

s 

E 

T 

HAPT, HPGT 

- 

Height at apogee and perigee for 
target orbit (km) ■ 

None 

T 

3 

RELINC 

1 

Same as set 2 

0 


T0, TT, 

TANOM0. TANOMT 


Same as set 2 

-1 


R0MAG, V0MAG, 
FLT0 

1%1- IVqI 
“0 

Magnitude of position and velocity 
vectors (km) and flight angle be- 
tween them for vehicle orbit 

ROMAG: -1 
FLT0: -1 

S 

£ 

RTMAG,VTMAG, 

l^tl l^tl 

Same for target orbit 

FLTT: -1 

T 

4 

FLTT .. 

“t 




RELINC 


Same as set 2 

0 


T0, TT, 
TANOM0JANOMT 


Same as set 2 

-1 
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C. Feedback Guidance Parameters 

In order to exercise the feedback guidance portion of the program and 
collect statistics on performance, the magnitude of the navigation update 
errors at the start of the first coast, at the start of the second coast and 
in the middle of the last burn need to be specified. The time between guid- 
ance updatefron coast and burn arcs needs to be specified and the number of 
separate Monte Carlo runs and time between statistical samples defined. The 
acceleration noise added at each guidance cycle is set at five percent of the 
thrust during burns and about 1/2 of the worst case gravity errors during coasts 
and can be changed if desired. 


Name 

Symbol 

Definition 

Default Value 

OELS(l) 



Time between guidance updates 
during burns (sec) 

20 sec 

DELS (2) 



Time between guidance updates 
during coasts (sec) 

100 sec 

NOISON 

• 


0 - no noise 

1 - navigation and accelera- 

tion perturbations 

0 

SIGI«R(i‘)TSIGMAV{ir 



Standard deviation of position 0 

and velocity navigation errors 
(km/sec2) at end of second from 
last burn (only used during 3- 
burn mission) 

SIGMAR(2),SIGMAV(2) 

'r- 


Same at end of next to last 
burn (km/sec^) 

0 

SIGMAR(3),SIGMAV(3) 

^R’ 

'v 

Same in the middle of last 
burn (km/sec^) 

0 

PERT(l) 

'a 


Standard deviation of accel- 
eration errors during burns 
(added each guidance cycle) 
(km/sec^) 

.05*T 

niQ 

PERT(2) 

6 


Standard deviation of accel- 

.5 10'^ 

a 


eration errors during coasts 
(km/sec2) 

R^ 

e 
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Name 

Symbol 

Definition 

Default Value 

MCARLO 

- 

Number of Monte Carlo cases 
to be run 

1 

PTB 

- 

Time between output samples 
during burns (sec) 

10 sec 

PTC 

- 

Time between output samples 
during coasts (sec) 

100 sec 

D. General Parameters 

Included here are the 
NAMELIST NAMLSl input. 

remainder of the parameters which 

may be set by 

Name 

Symbol 

Definition 

Default Value 

N0TAR6 


-1 Targetting only 

0 Targetting and feedback 
guidance 

1 Guidance only using para- 

' meters read in by NAMELIST 
NAMLS2 

0 

NAVOFF 

■ 

0 Convergence status printed 
whenever output sample is 

“taken in gu i dance"nibde“ 

1 No print 

1 

lOUTPT " 

- 

Integer parameter defining 
output device 

6 

EERROR 


If eccentricity less than 
EERROR it is set equal 'to 0 

.01 

TERROR 

't 

If tug or target within this 
time (sec) tolerance of node 
or some mean anomaly, consid- 
ered at node or mean anomaly 

10. sec 

RERROR 


Differences in angles (relative .5 deg''ees 
inclination, etc) less than this 
tolerance will be ignored 


OBLATE - Weighting factor used in setting 0.0 

oblateness effects (subroutine 
PERTO) 





r^ame Symbol Definition Default Value 


AXIS(I) 

Axis of rotation of the earth, 

■ 0 1 

must be set in relation to co- 

0 


ordinate system chosen by target- 
ting when oblateness is activated 

1 

E. NOTARG-1 , 

Guidance Only Parameters (^NAMSL2) 



The following parameters will define the orbits of the target and 
vehicle, initial mass of the vehicle (all other vehicle parameters are set 
by NAMSLl or default options), the initial costate vector and times array 
needed to define the burn and coast arcs. 


Name 

Symbol 

Definition 

Default Value 

NBURNS 

- 

Number of burn arcs 

2 

X0, T0 

Xq» tQ 

State of vehicle at start of 
mission. tQ time at start of 
mission. 

T0==-1 

XT. TT 

Xt» 

State of target at time tj 

TT=-1 

Q0 

ia 

Initial costate 

None 

AMj3 

m^ 

Initial mass 

28803.1155 kg 
(63500 lbs) 

TIMES 

- 

Array of times defining start 
and end of coast and burn arcs 

None 
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Program Output 

The exact program output varies with the setting of the output control 
parameters NAVOFF, PTB, and PTC. The nature of the output, by subroutine, is 
as follows: 

® MAIN - error messages, impulsive approximation summary, program 
notes of convergence status, and converged targetting summary 

° AUXOUT - summary of current convergence status 

° PHASE - error messages, orbit-type message (e.g. ''CIRCULAR/CIRCULAR 
INCLINED ORBITS'), coast messages (when states must be advanced until 
proper phasing exists), and phasing-orbit messages (including rela- 
tive geometry, "desired" phasing orbit, and allowable phasing orbit) 

° GLMNTS - orbital elements and designation as to whether they are re- 
presentative of state at start or end of a burn. 

® STATIS - Monte Carlo summary 

° USTAT - state, costate, and magnitude of costate vectors 

PTB and PTC control the sample collection times in guidance mode, and NAVOFF ■ 
controls the shutoff of the convergence-status summary (from AUXOUT) during 
guidance mode. In addition, there exists an internal program variable, LPRINT, 
which when set to 1 produces voluminous output on each call to GUIDE, detailing 
state-and eostate-a-t-predef-i-ned-t-ime-s-on- each-eGast-arc~and--orbi-tal-elements at 
the beginning and end of each burn. .Because it gives so much output, and is 
unlikely to be needed over an entire run, IPRINT must be set within the program. 



Interdependence of Subroutines 


Note that the dashed lines indicate further calls which are adequately des- 
cribed in the GUIDE document (except for the addition, in subroutine GUIDE, 
of calls to UCOAST and GLMNTS for output purposes). 
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Subroutine MAIN 

A. Purpose 

The MAIN routine controls the overall operation of the targetting and 
guidance program. It has four major sections. The input section, which reads 
the input data '^escribed in a previous section and calculates the orbital and 
vehicle paraneters needed to perform the targetting and guidance; the phasing 
and impulsive-initialization section which determines the number of burn arcs, 
rotates the target orbit into the vehicle orbit plane and calculates the planar 
impulsive solution for the orbit transfer; the convergence section which first 
converges from the planar impulsive solution to a finite burn solution and then 
repeatedly reconverges with the target orbit plane rotated in ten degree stepfe 
until the desired relative inclination is obtained; the feedback guidance sec- 
tion which exercises the GUIDE algorithm in a realtime guidance environment, 
continually reconverging in the presence of perturbations and collecting Monte 
Carlo statistics on the performance of the algorithm. 

B. Major Parameters (Input parameters discussed in Section 3) 

Name Symbol Definition 


HPG0,-HP6TvHPGX--- 

HAP0, HAPT, HAPX 

A0. AT, AX. AP 

E0, ET, EX, EP 
VAP0, VAPT, VAPX, VAPP 


Height" at perigee for vehicle, target and trans 
fer orbits (km. ) 


Height at apogee for vehicle, target and trans- 
fer orbits (km. ) 


a Semi-major axis for vehicle, target, transfer 

and phasing orbits (km.) 

e Eccentricity for respective orbits 

ly I Velocity magnitude at apogee for respective 
orbits {km/sec) 


VPG0, VPGT, VPGX, VPGP |Vp| 


Velocity magnitude at perigee for respective 
orbits (km/sec) 


TAU0. TAUT, TAUX, TAUP t 
IBOUND 


Period for respective orbits (sec.) 

0 - Outbound mission 

1 - Inbound mission 
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C. Method of Computation 

After reading the data (as previously discussed), the routine deter- 
mines whether a two burn mission will be sufficient. If the position and 
velocity vectors and the mean and true anomalies are not given, the true anom- 
alies are arbitrarily chosen such that a two burn mission is possible. This 
is accomplished by choosing the vehicle state, for T0=2OOO seconds, at node 
(perigee for outbound and apogee for inbound) in a coordinate system where 
perigee is in the x-j direction. This forces the first burn to be centered at 
2000 seconds and by choosing the target state at its opposite node (apogee for 
outbound and perigee for inbound) at TT=2000 + TAUX/2.0 (TAUX is period of de- 
sired transfer orbit) a two burn transfer is possible. For all other mission 

t 

definitions the PHASE routine is called and it determines whether two burns : 
will be sufficient and returns the state vectors defined at the time when the 
first burn is to begin. 

Impulsive Initialization 

An impulsive approximation is used as an initial guess for converging 
to the desired finite burn solutions. It is assumed that the optimal orbit 
transfer always has a burn centered about the greater apogee and this implies 
that the transfer orbit has as apogee the larger of the two apogees and as 
perigee the peri gee~of the other “Ofbit:; By calcuTating“the“VeTocities at 
apogee and perigee along the transfer orbit, the ay's required are easily de- 
termined. By converting these Av's to finite burn times, while assuming that 
the burns are centered at the respective nodes, and starting the mission 2000 
seconds before the node, a reasonable time history for a coast-burn-coast-burn 
mission Is defined. A reasonable estimate of initial costate is also needed 
in order to converge the GUIDE algorithm. By investigating the impulsive case, 
it is determined that the direction of thrust at the node is parallel to the 
velocity vector and that the rate of change of thrust direction is anti-parallel 
to the radius vector. (The reverse directions when decreasing velocity is 
required, on inbound missions.) By noting that the Iq^l is arbitrary for the 
boundary value, problem only one parameter. was left to be determined, the n- 
lationship between the |u| and |u|. (Note: qj = (u'*', u"*").) Using the fact 
that the variations in r,v form the same class of solutions as u,u, and applying 
the switching condition that |u| at perigee must equal the |u| . at apogee, it 
was found that the impulsive solution for u and u at apogee and perigee becomes 
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where r , v are the magnitudes of the position and velocity vectors at apogee; 
a a 

r , Vp are position and velocity magnitudes at perigee, and r, v are position 
and velocity magnitudes at either apogee or perigee (depending on where qg is 
desired) along the transfer orbit. In the program these formulas are further 
reduced and the |D| is chosen to be unit magnitude. The formulas become 



Cl = -fi . factor 


where for perigee the factor becomes 


FACTOR = 


0 * e/2 


e2/2). 


P P 


e “ transfer orbit 
eccentricity 


and at apogee it is 


FACTOR = ^ ^ P ^ 


When the mission is inbound and velocity needs to be reduced, the sign on 
both u and u is reversed. Since this qg is defined for the impulsive case it 
is good at the node and needs to be propagated back to T0, the chosen starting 
time for the mission. The two burn approximate solution is now completed and 
the program easily converges from this to the true solution. 

The approximate solution for the three burn mission is identical to 
that of the two burn one, except for insertion of a phasing orbit of pericd 
TAUP. For the approximate solution the phasing orbit is assumed to have the 
same perigee as the transfer orbit and the vehicle orbit (outbound) or target 
orbit (inbound).. This implies that, the burn at perigee is split into two 
burns and TAUP is chosen in subroutine PHASE to allow these burns to be aS 
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nearly equal as possible. The typical inbound mission approximate solution 
thus consists of an initial burn centered at apogee of the vehicle orbit, a 
coast from apogee to perigee along the transfer orbit, a second burn centered 
about perigee of the transfer orbit, a second coast of the orbital period 
(perigee to perigee) along the phasing orbit and a final burn centered again 
at perigee. The costate vector for the inbound 3 burn planar mission (plane 
change is added after initial convergence) was initialized using the samt for- 
mulas as the two burn case and the inbound mission successfully converges. 

For the three burn outbound mission, convergence proved to be more dif- 
ficult; It was discovered that the switching condition along the transfer orbit 
coast was very sensitive, and that the peaking characteristic of |u| at apogee 
and perigee was impossible to maintain when the phasing orbit was encountered | 
before the transfer orbit. It was found that by solving the mission backwards 
and integrating over the transfer orbit first, reasonable convergence was at- 
tained. In order to run the GUIDE algorithm baickwards from apogee on the tar- 
get orbit to perigee on the vehicle orbit with increasing mass it was necessary 
to make the orbits retrograde by changing the sign of their velocity vectors, 
to change the mass rate from positive to negative, to change the sign on initial 
u, to reduce initial mass and to alter the TIMES array. The TIMES array for the 
backwards three burn outbound mission 1s initially targetting by choosing it to 
be 

TIMES(l) = T0 

TIMES{2) = TIMES(l) + BURN3 

TIMES(3) = TIMES(2) + TAUX/2 - *BURN3 + BURN2)/2.0 

TIMES(4) = TIMES(3) + BURN2 

TIMES(5) = TIMES(4) + TAUP - (BURN2 + BURN1)/2.0 

TIMES(6) = TIMES(5) + BURNT 

Where BURN! is the length of the burn at perigee of the vehicle orbit, BURN2 
is the length of the burn at perigee of the phasing orbit, BURNS is the length 
of the burn at apogee of the transfer orbit and TAUX and TAUP are the periods 
of the transfer and phasing orbits respectively. The initial mass is reduced to 

. mg = mQ - iii(BURNl + BURN2 h BURNS) 

where m is positive. Q0 is initialized at apogee of the transfer orbit and then 
the last three components are changed in sign (u). 
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Mission Convergence 

Using these approximate solutions for the two and three burn missions, 
the planar missions are converged in less than twenty iterations. At this 
point the relative inclination, RELINC, between the target and vehicle orbits 
is tested and if it exceeds some minimum value, the mission is altered to in- 
clude the de'iired plane change. The target Oibit is rotated in maximum of 10 
steps from the vehicle orbital plane, and is reconverged at each step in the 
process. The two burn missions converged readily using this procedure but it 
was necessary to alter the three burn missions to two burn ones to obtain good 
convergence properties. This was accomplished by replacing the lowest orbit 
(target orbit for inbound and vehicle orbit for outbound) by the phasing orbit 
found during the planar mission convergence. The inbound mission is converge^ 
as a two burn one with the desired end conditions being the phasing orbit ro- 
tated about perigee. The outbound mission is converged backwards rotating at 
each step the target orbit as well as initial costate and converging to the 
phasing orbit. After inclusion of the total desired angular rotation, the 
third burn is again introduced into the mission definition and convergence for 
the three burn mission is attained. The outbound 3-burn mission is then turned 
around and solved in a forwards fashion using the final costate as initial co- 
state and the burn and coast times derived from the backwards convergence. 

Feedback Guidance 

At this point targetting is completed and a converged solution exists 
for guiding the vehicle into the target orbit. In the MAIN routine the major 
guidance function performed is to control the collection of and print the Monte 
Carlo statistics generated when doing feedback guidance. The routines BCBCB 
and CBCB called by MAIN add perturbations into the state of the vehicle and 
move step by step in time through a full feedback guidance cycle. At several 
points along each burn and coast arc, error statistics are gathered and an 
estimate is made of the error in meeting desired end conditions. These sta- 
tistics are collected over MCARLO separate orbital transfers and a summary 
printout is obtained from routine STATIS. . 

This completes the description of the MAIN routine. A math flowchart 
of it is contained on the next three pages. 



















RETURN TO RECONVERGE TOTAL MISSION 
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3-BURN OUTBOUND 
MISSION WHICH 
BEEN CONVERGED 
BACKWARDS, TURN , 
AROUND AND 
RECONVERGE 


IBACK=0 



RECONVERGE MISSION 


3-BURN OUTBOUND 

TURN AROUND MISSION TO 
INSURE CONVERGENCE 




REPEAT 

FOR 

MONTE 

CARLO 

RUNS 



ALL TARGEHIN 


260 ) 


G COMPLETED 


-Nk. 

PRINT PERTINENT 
DATA ABOUT ACHIEVED 
FINAL ORBIT 


NOTARG=-r? : 


GUIDANCE SECTION 



NO_ FEEDBACK 
GUIDANCE 
DESIRED 



STOP 


YES 


’^NO 


REMOVE ANY INITIAL COAST 
FROM 3-BURN MISSIONS 


NO 


IBOUND=l ?> 
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-G'CONVc-RT --HElGiiTb aT- -AND- PERI GciL --- 1 n TO -Ui*3-I TAL- LLE'*‘<Ei'JT.S»..— MAI 0-023-3 

bO 


3003 


AO-RuARTH-t { MAPU } /2 - 

EO- ( HAPO + KEAkI i-i ) / Al> — i .;■) 

• - 1 F- ( rj 0 • L T • »“ LR'RO W ) • L 0 • O — . 

AT — KLAi\Th+ ! HAP 1 +r’iPvj f ) /2 • 
LT = tHAPT+RLAi<TH)/AT-i .0 
— —- ••! F-(-LT'* LT *Ei_RRGF) ■ b-T— O . 0-- ■ 
GO TU UQ 


MAI00240 
MAI 0024 1 

MAI 00242 

MAI 0024O 
. MA100244 
— .MAI .0 0243 


C CONVERT POST r I ON* VLL.I )Ci TV ANl) FL IGH7 ANGLES 

.7-5 A0=aK4R0MAG/ ( ? * 4, i;<-VOS4Au4«24RbMAC ) : - 

H0MAG-DA2;> i R0MAG4 v ^rDS I N( FLY 0*1->EGUJN ) ) 

. . tO-OSOKT ( 1 . I A0>UX ) ) 

— .-.I ...j F-i Lo ♦ L T • c E rror ) • 3 u =o-* c — — — . - • — — — • 

ATsU^4RTMAG/ ( 2 *4Uk.-V f MAo4«2 4RTMAC, ) 

HTMAG = DAr3S(R7MAo=i-VT 1AG«-0SIN( FLTT^1 su.l<.;C0.'< > ) 


TO 0R3ITAL ELL'^EnTS* 


•VJAI 00246 
MA10024 / 

MAI 01x242 

MAI 00244 
MA I 0U23O 

M A 10-0.23; 

MAI002S2 
MAI 002S3 


LT=;DS0RT ( 1 .-H I MAG^'«::«-2/ ( AT^.SUN) ) . . . . - _.MA10U23^ 

IF ILT.LT .LE'rRDR ) LT-U.U MAI002-OD 

... 100 HAPu=A U^i' C 1 .+EU )-UEA;i r H MAI 002S6 

HPG0=A0^( 1 ••-f;0 ) -i<i.ARri{ ' - -- - •• ^ MAI,Ow2b7 

HAPT-=AT*( i -+E TJ-KL-.ARrp MA1U02S3 

h!PGT-ATJSt( I .-FT ) -rear r P MAlOU2Sv 

I F.( HPGO *GT , 0 . 0 • AnO. F.^’G r .gT *0 .0 ) . GU TO 110 - . .MA100260 

Wil ITE I lOUTf-J r *3i SO ) MPG0*5»PGT MAi002ul 

3150 FOkMaTC* cXECUIlUii i RAT 1 WG • • • / . ' Ht_lG.lT AT PLRIGLE UF TUG~* MA1002G2 

1,0 1^*0 * / • ■ PC 1 A r i^EKi-GcE -UF — TARGET- * »U 14.5) AI -0.0263- 

STfjP Ma100264 

C MAI oo 263 

C- UETL'.RM 1 NC INITIAL i.ND FIRAe- i^ADI !• TEST FIRST Ir AN iNuUUN MAi002oG 

C OUTEOUi-iD MlSSluiJ- li- APLAjLi'. uK THE rAKGLT QkfIT IS LLSS THiN AFUGEE .MA10J267 

C DF The tug UK.iiT, A.J IRGi-'Uriu MlG5iO('< IS ASSU:-‘3D. MAi002Gi.2 

C—CHECK-FLiR.-Gi<CAT;«R.../Fu.^L.x* — - MAI-0_02o3 

110 IF (HAPO . G r , I lAI’T I GU TU i2i> .MA10O27o 

C OUTCGUND MiSSJON (FROM P^-RIGLE TO APOGEE)* MAI0u27i 

IBUUND=C .. - . .. MA1002.72. 

..-.. RI=REARTH4HPG0 ' ‘ MA10U273 

rF = REAKTH + HAPT .MAI 00274 




MAIN - 14 


FI LtL I M A i 


r U ; % r ! ^ A ! 1 


CAi'^riKlDGL: ML)AiTO SYSTtM 


-O-lNBUUNU -W I ^>r.I urvi. if kUM -Al^uGLe-.7Ai -PuAiGLi- i .- 
I2b ' Ui{.UjAU=l 

Rl=.^EAj<TtJ + l lARO 

= ARTH+i iRGT- --• ■••— - 

C-CALCULATL. iNiTiAL ANL' h-IKAL VLLU C 1 Ti A:. S , 

.-ISO VI. = L)S(JR r ( (2 ./K 1-1 •/AG )-) - . - 

VF-=USQRT( (2 ./Kr--Tl.*/A.T.).) * - — 

C OCTERMINL VLLOCITI 1 .S AT ;.fwG POINTS Or •UGtTARGcT 
VAPO^OSOKT (OK-" ( ? . / ( iiAPu+RFAF^TH ) -1 ./AO ) ’ 

-VPOO^OSQRT ( UK^r. ( , /( 1 sJ-OO+KEAR TH )-l ./AO > ) • 


ORBITS 


•VAP'I = u SCJkT C ijK-!' C ;^ . / ( i iAP I + Tm > — l'« / A T ) ) - - 

VHGT -oSOk T ( UK v ( 0 . / t Mi 'Cl +i<LAR TH ) — i . /A 1 ) ) 

•O—OcF-I-NE- Si ATES-.iN-ARK 1 T.4 AA Y- -CUUKD-INaTlO-. SY OT 0.M *- UNLESS . UNE-. i- .IMPL^I-ED — — 
-C- THKOUGH- THE STATL VLCl'ORS ( PHASE V/AS OALLdDJ. TiiE TUG WILL BE LOCATED 
C AT PLR IGtE /APUGL.L And JHL TARGET AT APuGtE'/Pi-R lGt.F. DEPLNDKG UN THE 

-C— VALUE UF- r BOUND. P FR 1 GF.L 1 S - I N - THE DI RuC Ti UN (1.0.0) AND T E . H 1, 

-C— VECTOR- iNj. THE DlRuCTluN (0*0.1) 

X-l.F— IPHASL= 1 . THt- STAiL V,.CTuRS V/LRE St_T UP IN SUBKiJ’JTlNu PHASE. .. ... 
1 F-(-IH,HASE.EU-.1 ) . wQ -TO • 1 7u- - - . - • 


^5AIO-O27a- 

MAi 0O27T 
MAI U027B 

...MAIOGET^ 

MAI 0 OEOO- 

MAi.0C2E 1 

MA i.OOEBE- 

M A 1 0 0 2 B B 
MAI 0020^* 
.....MAl002Eb 

. -MAI 0 OEbO 

. MaI 0 u 2 S 7 , 
-MAj.aOEES 
MAI Q026^> 
MAI 00200 
_MAI 0 u.i 9 1 


S1GN=1 .0 
' IF ( IHUUND. LO . 1 > 
~KO (-1-) =Kl »i*S10N. - 


S I Gl'i— — 1 * G 


-C 
.175 


iTO(Ei-D.U 

RO (3J = 0.0 

VO ( IT-O.O - 

V0(2)=VI+SIGN 

vot:-iJ = o.u 

RTt l ) = -RF--J=SIGN 

RT(2)-0.0 
... . KTl o) = 0.U 

VT4 1 )-0*o 

vr(2)--VF:{»SiGN 
VT (3 )=0.0 

SET— UP '1 RANSFER ORril 


. .. A A — ( R 1 + K F ) / 2 # 

t.X = DMAXl CR i *RF ) /AX-i . U . 

— HAPX=AX;Ji(l .-fLX J-uLA.TTH. ^ 

HPGX = AX:i‘ ( I >-K'UAR TH - ' , 

VAPX=DSGRf(UX;j'(i’./D-<AAl(RI*RF>— 1-/Ax>) 

-VPGX=DSGRT ( UK* ( 2 . / D 'W r.4 1 ( R 1 «R F ) — 1 . / AX ) ) - _ 

-C-DEFiNE DELTA V»S (+ JR -) 

.-C-THE. INl U AL DELTA \/ iS Ui-FINLJ AS t HE .. 7LL0 C 1 TY AT Thi_ APuG X/PlK IGEE 

-C-ON— THL--TKANSFi_R ORE i, T . M i i ,US. -I H. INiTl-AL VELUCITY. THL F 1 NA.L . DELT/v_V._ 15HA1 0 Oo.l c? 

C DErliMLD AS THri FINAL Vi.LJClTY MINUS fiL: VELOCITY AT PEKlGE /APUGLE QN MA100319 


...MAI 00202 
, MAI002yj 

MAl-O-UEU.A. 

MAI 00295 
MAI 00296 

MAI,QCJ2y7., 

. ....jMAI 002.9O 

MAI00299 

MA.L003UU . 

MAIOCSuI 
MAI 00302 
MAI 0O303 

MAi OO lOa 

M A 1 0 0 J 0 b 

_MAi O.OEL 16 

MAI oOSo7 
M A I 0 0 3 0 6 
MAI0U30U 
MaICOoIU 
MAI 00311 

MAI 0 0.3J.2.. 

MA10o3l3 
MAI 003 1 A 
_ MAl 0031 S. 
MAI 00310 
MAI 0031 7 


C THF 


TU - 1 iiO 


transfer ELLlPSi:. 

-I F- ( lliUUND .LO • 1 ) -kiu 

Ut.LTVi=VPGX-7l - . 

DElTs/F =VT'-VaPX 

^VELXFR=VP0X — . .2- 

GO TU 19S 

160 L>LLTVi=VA!->X-Vi 

DI"LTVF=-VF“\/P0X - . 

VELXFR=VAPX 

,C._ OETEK^iNE. 0::El'kAL r’E;;ijDS* 

-195 TAU0=2 .*Pi *DS0R T ( AG**3/UK-)- 


MA100320 

_VA10032.1,. 

, MAI 00322 
MAIOO’323 

dAl .0 03 Z/L. 

MAI 0032*3 
MA I 00326 
M AI003?r 


. MAI 0 032 6 
-MAI0032V 
M.A 1-0.0330. 


MAIN - 15 


file: main 


f‘ iJi J Tl< Ai* 


1 


CAM..',h: lOGL MUNITU-^ 5YSTLM 


A — It:'-* 1 i ( A A '■? -i'-A/ iJ/wJ - 

TAUT = ? *^p I !^i>S(vR'n A ) ' 

OETEiihMUr. finite NUkNS FnU'.’ Dt! ?il PLO OF.lTA VS AM.) INITIAL NASS, 

Dunn l-=^A.’.U)/A.v.iJi.n ■*-( fJi.XFi — aMDO T *NALiL( ULL TVI J/THNUS T > ~L*.0.-i 

UllT AM -UUKN i :?AMDUl 

UUkNL=“( A.Mu -Lti LTAM > /AMUUT^ ( UE^XF I -AMLDT *L>AbS (DELTVF)/TH .UST)-1*0) 

CUALT-1.— TAUX/t: *-=r.Ci>U/<N.L+t_UicWL J v'a-* ^ 

SET UP TIMES ANHAY rUP A-UUi^N MI.SSINN WITH 
SECONDS uEFOFL Trit. *» UC IS UUE. AT Tf'Jt; i-iUDL# 

( IFHASL = 1 ' 

IT IS - GL . 


-MaT-w G^— • i 
MAI 

Mai 0033- 

.^3A>LC^.033S 

MAI 0Q33S 
MAI OOSoo 
J?lAijElO,33Z. 


C 
C 

•C-C3P IGINALLV SP tl C I F ,UE; AS- - GL . -U *0 
C RUN bAGK bY 5bC» 1 -k.» LONG AS 

T I M L S ( 1 ) = 0 . u 

IF( IPHASS.LGj.O ) Gl) iO IGo 

F r=o .0 

lF-( T0^i.'Gbb **LT*G*o> SMlFT=AuOU •~TU 

3050 FORMATE* TIME SCALL LrlANoLO bV.FlO.^** 

,..1AL COAST.**/.' TO ANO TT NOW EUUAL 0,0 

IF(-TO“3-0.uO-U*LT.C*U TT-TT>2UUC*-'TU 

1F( Tu-^000 *0 ,lt*g * 0 ) ro-0.0 


TIME TO-U ARblTfTlARlLY 
UNLESS TO AND T T WERE 
IN 'a'HLLH CASL . . iO„.iS-.. 

O • V « 


2000 MA100333 
MA10U339 

MAX-UQ340 

MAI 00341 
MAl 00343 
M A 1 Q 03-4^ 


slcoinus to 

AfiD * .F 1 0. 2 ) 


ALLOW FUR 


196 


IFISMl FT ,GT . .00 1 ) vN-v i TE I I OUT P7 ♦ 30 50 ) 

..•IF.ITO.GT.0.001 ) -TO- J 0-200.0* 

GQ TU 1V7 
T T “20 0 0 • t- T AUX/2 • 


SHIFT. TT 


M A I 0 0344 
MAI 00343 

MA1.QU340 

i N 1 T 1 M A I 0 0 3 4 / 

. , MAI0U345 

L‘lAl-0.0340 

MAI 0033C 
-MAI 0035i 

,MA.l.U.Q3.52 

, MA1U03S3 
MAI 003S4 
MJU..0 a.35L 


197 TIMLS(3)=if0c0.-NURNi/2.+T0 

T IMESE 4 )=1 I MuS (3 ) +r;URM 

TlMtS(.5)=TIMCSI4 ) +LLIAST.I: - 

. TlML5tO)=TlMLS(S}+f‘Uf:N*_ 

C.. STUG-iu THE TUG ANu J'.ARUbT SfATbS 

UO--2CO- .1=1 »3-- 

,X0 ( 1 ) =f<U( I ) 

X0( J+3)=V0C i i 

XT ( I ) =RT ( 1 ) 

200 XT(I+3)-vTU) 

.C.. DEFINE IMPULSiVL 


IN r:lL WORDING VARIALLLS X .Xf 


MAlOO^So 
MAI 00357 
._MA 10.0355 
MAI 0035*^ 
MAI 00360 
_-MA10036 1 
MAI 00362 

MAI003o 3 
_.MAI U0364 
MAI 0036s 
MAiOUSoo 
_MAI.aO„3.6.7 
W A 1 0036 5 
MAI OQ^oU 
„MAI0037C 
MAI 0037 I 
MAI 00372 


CuSTAiL. U IS OF unit MAGNITUDE AfTO ALONG VLbDCiTY 

.C— VLCTUR .-A.i-'-.THL.-NuLiL.. -ANG... U -U J-T.-...Lu„.ALUNG . ENL-LAu IH RAO 1 US-.VLC. JJR* 

IF ( liinUND.LU* 1 ) GO TO 2C0-.> 

SlGN=l-0 

FACTUR = ( 1 *+-i-X/2 . -EX-:«-;^2/2 . ) -^UK./ t P i ’J^-XJ^VELXFR i . 

GO TU 2007 

2005 SloN=-i.O 

HALTE.)R = (UR+-\/AP.X^^..VP.oX.'r; I HAi-'X+REAkT H ) ) / ( ( HAPX+iiLART H ) '-^>3 -7^ VAP_X+_VtiGX.)J:MA1.0a373 

2007 CUNTlNUE MAIU037^ 

DO 2ui 1=1.3 MA10037S 

OOl.l)=XO(l +3 J / V i •i' S 1 GN - . : ... MAI 0.037 o 

201 00 ( I +3 ) =-X O' I I ) '^r AL VNI ;N.S 1 oN MAIC037 7 

NU = 0 - MAI00J7L 

.C-CaAbT--TUG.-3ACX...AR3 I li^ARY *_.0 00_...3L G. . ,uE.FURE...ST ART...uF . , rHL„.i„ST„.b.Ui<N* ^A.1.0,03'/:U 


* . Xt 


CALL COAST ( Xo • GO . -2o00. ,XO *O0 * PHl .PHI ) 
T3LG IN-T IMES (3 ) 

..Ga-TU...202 . - .. ' 

'X 'i^ T hi-\L L“L>Ui<N i.iiSSlliN ^ 4= 7 X*. . 


i=-ijURN- -'AiSSi 3N . LS- ASSU.'-iuU ...-UNLs-SS- ASi.NG -WIl-L- AeLUW. 


MAI 00330 
MA I 0038 i 

M A I 0 0.36 2 

M A I 0 0 3 3 3 
- M A 1 O 0 3 if. A 
,__MAi.0.u,3G b 


MurUTO :iYSTt_M 


— C — A -- id— ,-M I i>S I Ui» • - I r .. fi ii_ ...TU\j .. /'.i'» l-> T.AK Uw T AT-L.lj, . W wK£d...i.-P.£i.C.iE-it— D • — Ai. .0 u-i •— o 

Soo NbUKNS = 3 MAlG03t>7 

IF { POMAG T .0 . u) GO TO 3ui MAI0036 

PCWAG^OSOk'T ( AO ( 1 ) j;- 5-*i^+ .iO ( 2 )*S==X2 + A0 ( 3) > - — - JAAlii.o3as?. 

S/OMAG-OSOK r ( 0 C 1 > ❖•i‘3+\/0 t 2. j 4-4s2 + V0('3) ) - — -- - ,M A1 0 0390 

RTMAG=OSURT (AT ( 1 ) + K I.(2J**2+R.T (3)**2.) ... - - — MAi0 039 1 

^-V-TMAG.=DSOrM.C.V-T..( I ) ( 2.) >,'^>^2+VT-(..^>,.:4v5:2.) WAJ_aa393 

IF (ROMaG.G r . Ah A,‘-<TH. Ai>iD, VOV'AG *07 • . 1 . AND . MAl003h*3 

I RTiViAG*GT . ftLiA=-^TH. Af*L>. VTMAG. GT . . 1 ) G O T 0 6 0 0 MAI 00392 

-C-STATC5 -W<_Rri MLH ..GPcClFlGG* CH£C*< IF URU ITAl^ . LLLMGNTS WLRB. -„.MAl0039G 

ir ( AU • GL * kGAI^ITH* A03 . 2 O . GL • O . U *Ai-iL> • MAl 0039 6 

1 .... AT. Gi;.i>;6AATH . AGO •.._!. 06 . 0. 0 ) GO T(.i 310 .... „... . ..... .M AX 00^V>‘A 

-C-CrtECK~I-F~lii-:i-GHTG. A T -APUuGt- ANO - P-6R 1 uGG - I K.’4 ) -k/LRC.- SPt_CIF.lGD, 6d;U-G0-196 

I F ( HAPO.GL .0 *0 . A is O . 3i.»GU . oL . 0 . 0 . AND . H AP T , G3 • 0 . O MA I 0 0399 

1 . AND. f7PGT .G6 . 0 . O ) Gu T U 6o0 MAI 00206 

CHr. CK- IF- - MAGN I TGO^S DF_P;iSI T ION . AND VGLUCI T V ..ArjD -FGIGH L' MALOo.201 

-C-ANGGES. SPGCIriLU.C i- 6 i Oi 3 T . A NugL UGF 1Nl_D IN 2-uUi^N CUM'^uNTS*) MAIQU202 

-501 IF (POM AG .GF .KL- A?< 1 H. AND. VOmAC .uG. .1 .AND . FLTO . CG . 0 . 0 . AND , ...-MA1 00.2 0 J 

i-RTMAG- oG-.Kc A4^ 3' H . AND .. V-T .'0 AG *G2 . .1- ..A ND. F L-T T . GG - 0 . 0 ).-,GU. --TD,50:7 MA1 jJQ2 02 

C INADGOUArt GLGMr.i^TG GPGGlrluD. ST03->. MAi0C2O2 

WPi fF( I0UTPT»3 ioO) NO . \/0 .PT * Vr .RGLINC. AO .60 ♦ AT *c T MA1002Cb 

-31GO---FOPMAH 6 tXr:GUTJUN TGPM i N A T-i I'iG . I fflPPUPGP . DATA SP EC IF i iDO -FU P, . 3 -U.UtiNMAl.O 02 0 / 

; L...M1G31U3>*. • ./. * GwT 1 Ro = * . 3D i 2 . O . • V C- * . 30 1 2 . o . / * ’ kT = • . 30 1 MA 1 0 02 0 3 

S2.G»' VT-* .301 4. G./. • iiLulNC-* . L'12 .6 .// * * SET id AO =. • • Ul4 . 6 ■,MA i 002 09 

3-»-c0=*-»0l2.6.-*.. aT-» LI .014.6*/ ) MALOOALO 

GU TO 26 MAI 0041 I 

C MA100412 

-G-CONVCRT-HF: iGJ-iTM • iNTi.i •UiiG 1 T AL • LLLMG3-ITS . - - ...... ... MAIQO 2 I 2 

-bOb... A0 = l^LAPTr1+ U-JAPu+HPGO ) /2 . 0 ..MAl0u4l2 

CU-(HAl^O + PLAPTH)/AO-l .0 , . MA10041L 

lF-( LO.LT.LriAKQlP )• t lo-u.O- - J ■ .MA1,UU4 1 L 

AfaRLAPTH+ (HAP T + RPG T) /2 « o' MAI 004 I? 

hT- (HAPT + K'cAKTH) /A r-l . 0 MAI 0041O 

..lFU:T.LT*r;L-.Pi;03O LT = 0 .u . . . -. ... MA100419 

..... GO TO 610 MA100420 

..C -CUNVLl^T. POGi T 1 uN * VLLuC 1 T Y ANu FLIGHT Y.MGLLS iNfO OPLITAL LlLMLNTS., MAX0022 1 

-30-7- AO-UK^}<ONiAk>/ { 2 OMA vL P OMAG ) — MALO 0.2.62: 

H0MAG=DA)-i3( ivOi'-l AG-^Vv‘i'.i,Mj;>D5i i ti C FL T 0 ❖ DLG CQ3v )) MAI00423 

LO = DGOKT 1 I . — hlOMAG>2-.! 0/ t A 0 J ) MAI00222 

IF ( GO. LT .Lt.JGTOi-< ) t-.0. = o.O - MAl0022o 

AT-ON’f'PT MAG/l 3 f I'MAG ) MAi0022(j 

HTM AG— DAijG ( PTM.-\L» -1= VT i-1.\ lo 4 OG 1 Is I i" i_ T T * DLv> CO )) MA10022'/ 

E.T=D6GPT (-1 .-HTMaG^42/-(.A]>LI,U ) MA10Q4.L9 

iFIET.LT.LhKSvOP) FT=0.O MA100429 

510 MPGO = A 0^ C i • -LO ) -KhAA' m MAI0043O 

- -HPG T-A T=Jt ( I . -J.T ) -r,‘L6p T3i - - - -- - MAl0u2c.l 

-C . .• . .... MAI 00232 

..C...ELE.MLN7 S i'i LRL..., SPLC 1 F 1 Ld . Ci i^G.n . MLTNLP .TPOtl ANDMAUILS WLRL SPLCIFiLD. MAlO043j. 


-G-CONVCRT-HF: IG^■iT^.i • iNTi.i -UiiiG I T AL • LLLMG3-ITS . - - . . ...... ...... 

-bOb ... A0 = l^LAPTr1+ U-JAPu+HPGO ) /2 . 0 

CU-(HAl^O + PLAPTH)/AO-l .0 , ■ . 

■ ■ ' 1 F* ( LO • L T • L r;*< KOIP ) • i. o — u • O — - ~ ■ . -- . ■ -... ....■ . — - — 

Af=sRLAi<TH+(HAPT + KP )/2*U 

ft- (hapt+plaktm) /A r-i • 0 

..lFU:T.LT*FL-.Pi503O LT = 0.u ... .... .. 

..... GO TO 610 

..C -CUNVLl^T. POGi T 1 uN * VLLuC 1 T Y ANl^ FlIGMT Y.MGLLS iNfO OPLITAL LlLMLNTS. 

-bO-T- A 0— UK^}<OM Ak> / 1 c • 2GN ~ V O M A »L* P O.M As> ) — 

H0MAG=DA)-i3( ivOi'-l AG-^Vv‘i'.i,Mj;>D5i i ti C FL T 0 ❖ DLG CQlv ) ) 

LO = D60KT 1 I • — hlOM/\G>A-.! 0/ t A 0 J ) 

IF ( CO. LT .Lt..PPOi-< ) t-;0. = o.O - 

AT-ON’f'PTMAG/l 2 .7G,'s-V f :-^A.j-?72^i<rMAG ) 

..... HTM AG — D AuG ( PTM.'\L» -1= VT i'1.\ 4 OG 1 N I i" i_ T T ^ L'Cv> CO ) ) 

fcCT=D6GPT (-1 --HTMaG^^C/^.AT^LIP.) ) 

iFIET.LT.ChKSvOP) ET=0.O 
510 MPGO = AO^ C i »-LU )-KhAG' m 

- -HPG T— A T^ ( I . — t ) — r^i:. a; s { H - — - •— -- — — • - - 


..C...ELE.MLN7 S i'i LRL..., SPCC 1 F 1 Lg . Ci i^C.n . MLTNLP .TPOtl ANDMAUILS WLRL SPLCIFiCD. .MAlO043j. 

IF-CNPGg.Gg . 0 .0 - AND .Ki^G'l .-gC-O ..0 ) -..GO— 7 0-61.1 ..^AJ.U-0.23-2. 

WR 1 TL ( iOUTP r , 3 IGO ) HPGO.HPGT MAIOOAjS 

STOP MAI004GO 

6il.l IF (-TANOMO -GL-0 • O' . AtsO . TAND.M T- GL.-O .-0 ) GG-TO-660 MAIO.C.23,7. 

X~-ANQ(MAL.1L6.. NOT 6P L C 1 r 1 C D « . . T hi L TIMLG T0»TT vMLL f«OW Bu UN3 IDLHLD .. MA10023-3 

-C — ^THL- MCAW A Nil‘*‘ A C I t_G ( 3 iI*UG iJiNCA. PCNlCCu. OP* li* CiPCUCAP* G .NCL. T HL -MAi0u239 
-C—NOD'L) - ' - - - — — - 3V3Al.00.22o 


MAIN - 17 


FI LE I MAIN 


r OK 1 Is! Ai-i 


1 


LAM;;r-<iOGu MaNirQ.-. SYSTEM 


Ki^l-TAE 


1 »T-. E-i Tr-JL K - -1 i. • i>K w f i . F i iAi'i— -1 I ■ 

T AU 0 = i? • I ^OSOhT { ) 

TAUT=Z .=s^p I ^OSUKT { A \ ''r^ '^/UK > 

1 F-(.-TO * L r • TAUO • ANi> • TT . LT—T AUT- ) -Gn-tO--52 0 

• C TIMES AHii NOV LESS I HAN iJtMu UkUlTAL PLPiOD* STOP 

- WK1TL( 1l}U 1PT*J10L ) 7 0 ♦ r T , I AU 0 • T AOT 

3H>2 format (-* —LaLLU T l ON- T LRM 1 Na T i NO AN . ANOMAl. I ES- EXCEED-. 

I PERIOD**/** r U= • . L> i 4 . o , » TT-»,Dl4*6** T AUO= * * D 1 a • E * » 


>I-U.L 


16) 

C TlMi.S SINCE PAoSAGl; mRL REASGNAiiLE* LErlfJE 

C-STATlS at HEi^iCLL A.'.D HjNjPACATE AHEAU VIA CALLS TO 

■ C ~ ■ T -TME~S-~ C I 'VL N - t> Y F HE- .-AN - ANtiMAL-I L S * ' ■ — 

520 STATE! 1 >=A0^ ( 1 .-EO ) 

STATfc(2)-0,0 

STA'fr(S>-0.0 - - — 

STATL<E)=u.U 

STATE { S) =DSC»U (UNA I 2 . /STATE ( 1) -1 . /AO ) ) 

N0=-1 

CALL COAST(STATl. ,0;)UM* iO*srATL,OOOM*Phl*Phn ) , 

--00- 626. I = i *3 - 

fiOt I ) = STArE( u 
VOl D-STATEti + J) 

—STATE! I-)~A r^( I .-uT ) 

STArE(2)-0,0 
STATE(3)=0-0 

— VMACNTaDSOR » (UN* ( 2 ../ >Ta TE ( 1) -1 ./AT) ) - 

STa7E(4)—0«o. 

S T A r L '( ) — V .’"i A C N I O ^ S ( K I'i L 1 N’ C D E C U N ) 

—ST A T E ( O. ) = V A ON T L S i v * ( i 4 4 ; c 1 NC-i^ 6 E L C □ N J_ 


525 


MAl0o442 
MA100443 

y.A-l-OO.4^4.. 

MA100446 
. . MAI0044O 

-NL— UR-6I.TAI JAAJ-Oo.4-47. 

:AUT = * *D14, MA10044S 

MAI00449 

-MAI 0L4::iS 

TUG AND Target maioo4Si 

COAST* Q THE MA100462 

MAJ, U 045S 

MAI 00454 
MAI004bb 

MAI 00456 

- MAi0045'T 

MAI O04bo 

«__MA I 00.450 
iMAl 00460 
MAi0046l 
-_MA 1-0,0462 
,, MA100463 

.MAI 00464 
..MAX 004uS 
MAI 00460. 
MAI004o7 
.-MA1.004CiS 
MAI 0046/ 
MAI 00470 
-,.>^Ai.004.71 


*06UM*PrtI ,PHI ) 


530 
C- SET 

-C- 


states through T!LIR 


CALL CUAS r ! STATir. * 06U M * TT, STA T{ 

DO 630 1=1*3 

HT( I J=STATL( n 

VT( I )=STArEt 1 + 0) 

bOTH TlMuS Tu 200u. ArTE.'< PRUP/^GATING 

ANU.MALiE;»- - ^-v-: 

T0=2000* 

Tr=20D0.0 

— GiJ -Tt; 600 - 

_..C- 

.....c. TRUE ANDMALILS .vERE Sh.. e i r i El) • SlT UP COORDINATE SYSTEM SU 

C THAT THE - X!!) AXlS i S TU AnDS TU 6- Pi_R iGc.L • X ( S ).--I.S . ALONG 

C H* AND X!2) IS X(3) CUuS.'- X!l)* 

R«AG=A I 1 . -r.0=>= ) / ( i . + Ku*DCnS ( 1 ANuM O^DEGCON ) ) 

-RO ( 1 > =KMAG'-%UCUS { T ANU'T 0AiUv.-.GCON ) - - 

kO ( 2 ) — KM AG ^D S I N ( T A ^'4ij. "i 04 DuGCO N ) 

K0!3)=0.0 

— VMAG=DSURT (UN/ ( a 04 ! i.,.-uU442 ).) ) •• 

VO ! X ) ^—VM AG4DS 1 i-i ( T AImi iMu 4DEGC Ui'.j ) 

VO ! 2 ) =VMAG4 { EO +UCO'., ( TANQM0 4D EGCUN ) > 

.-VO!G) = 0,U ■ -. ■ - 

RMAG=AT4! 1 • -LT4 4;^ ) / ( 1 . + E T 4UC0S ! T ANDM T+ DEGCON ) } 

. RT ! i ) =RMAG4DCUS ( TAi;U‘l r-vDLGCON ) 

•-R-T- ( 2 ) =KMAG4i.»Si ( T A( A.ji-;T ^^Dl-oCU-N ) • - - - - — 


H 


550 


MAI 00472 
MA I 00476 
. MAX0Q474 
MAi 00475 
MA I 0047o 
„HA I.O 0 4 7 7 
MAI U04 76 
M.A10047/ 
.MAi.O 0460 
MAI 0046 1 
MA100466 

MAi.00.4.6.6 

MA I 00464 
MAI 004to5 

MAI 0048 6 

MAI 0046 / 
MAI 00436 

MAI -0.0 AES 

MAI 00*+9^ 
MA100491 
MAI00492 

MAi004/j, 

-MAIC049'^ 

- MAlOuASE 



fill:: mak; 


' clr-v T;'v A f * 



CAMiJKiaGL: MONITO-: SYGTLM 


G)=U*0 • \i 

VMAO — DSCi^if (OK./(/'»T-vl 1 m } J) 

VT( 1 ) = -vmag->-l>g 1 f i rrAr4r.-viT'-^L>LGCUN } 

VT ( 2 )-=.VMAGv ( L C +uLuS L I A^iU-.‘. T LGLUN J i 

VT(G)=C.O . -- 

-C— ADD—lr<CLlKAT 1 o.'-i TU -T AU^LT. aUulT. ... 

t^T(4 )=ii-T-( l ) V- 

KT ( 3 )=kT ( 2 ) ’!‘L>F IN { HiLL I NC^L>i£GCON ) 

RT ( 2 ) =RT< 2 ) ->:OCO‘: { Kl'L I wC'=J'Ut.GCON > 
VT( 1 )=VT4 1 ) 

VT ( 3 )=VT ( 2 ) vDS IN ( r':LL I ) 

. VT ( 2 ) =VT ( 2 ) L i<LL 1 »nL UN ) 


— T 0—2 0 0 0-* O 
TT s=2UG 0 • 0 
TD" 


MAIN - 18 


MAlOOAvT 
r4Ai0 04y- 
.MAi0 042v 
MAI OGbOu 
-MAIUObU i 

_MAmoi>a2 
MAlOObOi 
MAI OObOA 
...MAlOQbOU 
.MAXOObO^-' 
...MAI O0307' 
_JyiA.LUO.30i. 


•.v'LR! 


'LClFiLU. CHLCK iF THL TIMLS A.L T»C 
AT Tt-ii_ SAMt TiML ‘Vi-lLN PHASE I i.- .CALLED. 
CLUCN -T I.v.L, AND....NEVLR DL CRE ASlE-S*- 
iU 70 U 

G,< HUK.aARDS IN TIME AS RLClUlR-O. 
-UECJUR . Ur STATES* 


-G'O -TD" 700 

,.C-.-TUG AND TAKGLT states 
_C...THE- STATi:.S MU-ST iJh SPuClFlED 
_C-T-0— A-S~-ASSUME-D - TU THii-.Ri-AL- 

600 I F ( DAB S( ro — fT)*l-7».i) GU 
C COAST TARGET STATii ; i»-sCN '■^ARDS 
-C-SUaROUT INL COAST EXPECT;.* • A - 6* 

-. .DO. 6U I i = l *3 • • 

STATe( U=KTt I) - • 

-60*1 STATE t-I+3 J =^VTt I > ■ - — 

N0=-1 

CALL COAST ( S f ATI. , 0;.UJM * lO-VT . STATE * ODUM , PHI *PH1 ) 

DO 6UE l-i *3 

!U t 1 ) = .STATE( i) 

..602 V-T ( 1 J-STATE { 1+3 ) - :••— - - ; - 

C ARBIT.RAKILY CnANGL T i v,l ( JH 1 1,» I N , I F Tu*LT* 2U0U*» INCREASING 
C SECONDS TO AuuUW rUR HLiSSInLL 1 N I T I Ai_ C OASTS 

IF(T(>.Gr.2000. ) GO Tsi VUv- - 

TEMP-2000. -To 

: tt-tt + temp ■■■.■. 

WRlTuI IOUTPT,3120) 

3120 FUF^MAT(» CHANi,Li> Ti 

L*- INCREASED T k> AND 

_ 2E10.3* * * * »Llr>..i ) 

...CCALL PHASE TO DeTEk:-UN_ Ir UENU-.ZVOUS IS PUSSISLE .VlTH 2 UR 3 
aORNS»- FLirs:- C iRCULAN - T u AR - Cv'*"PL ANA P M 1 SSI UNS * A 2 —DUNN - - 

C MISSION, IS always F^ iSS i • FUR OHiER Gl^OMETPiuS* A 2 OR 3 BURN 
C MISSION may Bl. PO.SSlHLi-:* EXECllTION WILL T*iRa‘INATE IN PHASE. IF 

-C-MISvSlON IS - PUSSiSLL .7 1 f 1 1 1 N ..THE ■ ALLOTTED TIME. -- - 

__7q0 CAl,L phase ( :< 0 * VO »R 1 * V V . i'ibURwS * T CO as I *TAUP) 

iPHASE = l .... 

.__T: AUU.=2 *>i^P i 'SDSPf- l'.( Aw-1^-^^S/.L>n.) — . 


to TO 2UoO» 


U-Mp , TO. T T 
'-•ii., t SK i G 1 ru 
ri i.*Y **F-/.2* 


ALLOW FOR INITIAL 
• :.* TO AND TT NOW 


COASr.S. • ,/ . 
EUUAL .• .. - 


NO 


TAUT=2 *4;P1 >1^.030^ Tf AT ) 

Tu=TU+TCOAST 

_,.j:pJNaURNS.LD.2 ) CD ro 100 ..• 

-.c - DE TLU-M i NL. - V.'rlET H.-- .< AN ii'-.LA.>Ur*D OR OUroOUND MiSSlUl' 
— PORM-I NG" HE I CH TS -AT- A P t- A iNiD - PER I OllL-*- 


Is »<.EtiU 1 R lii . AF.Tf.R 


MAI 00300 
MAI 0031 0 

MAi 0031. j 

SAME.MAl 0031 1 
MAI OObl s 
:4Ai,0Q3L-' 
MAI 0051 3 
MAI 0051 w 

MAX 0.031 7 

»MAI Q 0 3 1 o 

.....MAI 0031 

MAI 0032 L' 

MAI Lj032 I 
MaI 00325 
_...MAI 00323 
. MAI 00324 
MAI0CS23 
._MaI 00323 
MAI 0052 y 
MAl 0052i' 

. - - M Ai O0520 
MAl C0330 
MAI00531 

MAi-0 us3s 

MAI C053.:. 
MAI 00534 
MAI0S35 5 
MAl 0053c.' 
MAI00S37 

^4AI.aS33E■ 

MAI 00550 
M A 1 0 0 34 C 

M A1.0034 i 

. -MA1U0342 
MAI 00343 

.._MAX.0.U.344 

MAI 00343 
MAI00543 

MAXUO.54 7 

. ...MAX0034O 
- MAI 0 0d 4V 
MA-IOOssu 


MAIN - 19 


FlLc; MAIN 


FOt'. f U Ai'J i ‘ t 


CA;.^OR100t£ (•lUi'inO/^ SYSTEM 


-4-»APv-=AO t +r.0 >-“.-^ uAN>ii 

HPoU=AO=i' ( 1 •-L-:u )~KilAi^TH 

HAPT = AT*( I.+LT H 

-HPGT-AT*( 1 .-clT )-;<i'AU i'H 

IF(HAPQ.GT*k-4AP1 ) • ISL!UIn1U=1 

AP= ( TAUP^»2^UK./ I 4 . I *3: 

- 1 f>U UND-- 1- 0-. 0 )_ H P ;,i • -H io-0 

iF( iPOUfND.EO. 1 ) hp<,h=hpgt 
EP=i .0~IHPCP-H’T:A>'TM )/Ai-» 


333 -.^3 


j-jAPP=aP*( I *+r.P - - - 

CALCULATE VLLUClTlES A> LNl> HUH-^TS Uf ALL GRl..1TS 

VAPO-USUPT ( U.<-f ( I • -uO ) / { ( 1 LO ) J ) 

! VPOU=USCU.T lUK;=l-(.i .+L-U ) /.(.A0;> C-l-^-LU ). ) ).-- 

V-APT— USUPT ( UK^ ( I • -L r ) / { AT (I . +I-T } > > 

VPOT=DSOkT t UK^ ( I - +1--T ) / ( AT^;^ ( I *-L T ) ) ) 

^YAPP— DSUU r ( UK'V: ( 1 , -i.P > / I AP« ( l-. + LP ) ) ) 

VPGP^USUKT ( ( i . / t AP-'-^ I 1 .-LP) ) ) 

IF( IBOUNO.LU.l ) CU 'lU yOl 

fii=hPGQ + «LAiiTl"i - — 


-7ai- 


.702 


RF=HAPT+RcA,iTH 

GO TfJ 70^1 

•-R 1 s'HAi-^ 0 +PE' Al?TM — •• 

KF-riPGT+KLAKTH 
V I = US0RT C U^'?‘ (2 */ii 1 . /AO M 

,.VF-— DSOWT CUK.J>. ( 2..*/r<K.-l.-/A.T..)-) — 

AX-CWl +«F )/2* 
tLX=DMAXt (P I .PF )/AX-i *0 

^.HAPX-AX*( 1 •- + LX ■ . .--- ■ — 

, MPGX^AX’S'l i .-LX )-i<LPl<TH 
VAPX-USOH r ( UK'S' ( ^ . /UMAX I (i< I tK F ) - 1 • / AX ) ) 
—VPGX— USOS5 I ( UK’^ CL' . / UK. i N I (*.U - i * / AX ) ) — 

IF ( I bound .c-G , 1 ) GU Tu 7o3 
DELTVl j:VPGP“Yi 


3dAI3>-0St>3., 
MAI OObo?-' 
MA I L»ubb3 
JAAI.G.OlSA- 
MA I 00 

MAlOObbO 


MA 1 OObbb 
MAI OObbS# 
_MAiO QS60. 
MAi OObu 1 
MAI 0056L 
__MAJ_O.U SB 3 
MA 1 0Qb64 
MA I Oubob 
.„MAl0ab6t_ 
... MAi 0 Oba 7 

MA I UObbS 

:AAl.o-uSbu, 

MAIOGb70 
MAI 00S7 I 
„><Ai.oab7,:L 
, M A 1 p Gb73 
,.MArOOU74 
—MALO U3.7bu 


.XP t -/.MDD T^DAUS (UuLTVl )/ IHPUSI ) -1 - O ) 


DELTV.M^VPGX-VPGP... 

. DtiLTVi-=Vr-VAPX 

GO TO 703S 

.703 ULLT VI- VAPX-tV i 

DLLTVM=VPOP-VPv».\ 

OcLT vf-vf -vpgp 

-C- F-l NO F I N I'i' L - k '.U/nN T 1 i . . 

703b UUkN 1 = “AMO/ AMOl.i r t u, 

DLLT am = a .'to U T =* rJ U.< N I 

bURN2=--( AMO-DLLT AM) / AM^O f i D LXP {.-AMD 0 T>sc Q AOS (OLL.T VMI/T.i-tj%US.T.J- 

ULLTAM=AMD«.JTsJ<bUR N2 + i>!.L T AM 

UURN3- — t AMO— Oi. LT AM ) /aW!>OT I DEXP ( — Ay4>UT *DAtjS iDELTvF )/TH *UST )— i.O) 

: IF l I 30UN0 .EU ..1 ) GO TO- / 041. - - 

,^UAST1-1 AUi^-(EUKNi+uUKS)S >/2. 

COASTk:=T AOX/C .-i COKn^4bUi<lvi3) /2 . 

704 COAST i -TAUX/E’ .-( OUkM + i .-U[PNP ) /2. . ' ■ 

C0AST2-TAUM-(ilUKNE + i AUvAO ) /2. 


MAl00b76 
MAI 00S77 

MAX.0,,Q5.7Li 

,MAlGub7y 
MAI OObdO 

MAi.QjjSUl 

MAlOObfc? 
MAI 00bs3 

MAiOuSciA 

MAi OC see 
MAI 0Ub3o 
MAI-OObb?. 


-l.OJ 


MAI OUbtio 
MAlOObbO 

MAi u.E.b^G 

. MAlGUb9 i 
MA100G9E 

MAi.OG.P.yw 

MAI 00b9s 
MAI U0b9b 

MAi0.1^3^.0C 

. MAi U0;>97 
MA 10 039.';; 

MAt.o.a.s.yE 

MAI 0 060 0 
MAI QOoO 1 


C SE T Ur:> -T iMLS .ARRAY - 

70 5 .. IF t TO -2.0 0 0 . .LT . E - 1; ) TT-TT+20 GO 

• iPlTO— 4-000. O.LT*0«u^ TO — 0*0 

j p { .T G • GT • O • O 0 1- ) - 1- O — T -O — 0 0 0 •— — — -.• 


-To 


MAI.U 0.6.03 

M A 1 0 0 u 0 j> 

, . MAI UOoOA 
MAX.UOGi:?^. 
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file:: majn 


F L>rv r F Ais 


CAMLiRlD&u MONna^ SYSTEM 


.-T0-=-2.U0U.-i>URN I/:. .+TO 

TiHciO I iN =T 0 : 

T1ML-S( 1 )= TO 

-TIf/T:^(.2 )-TO+i3URM • - 

TIM£2.( a ) = TiM:_2 C .c ) -I-i-Ua:*! 1. 
..T1MLS14 ) = T JMi_S t ^ ) +LiURN? 
-tl-MiiSt b ) ^T-l '4i:,vb L M > -huUA:;i-7 Z 


TI M£S( o ) = T iMLS ( ) -riiUn I'iO 

STORtr STATtfS IN AO Ai^jL; XT 

00-710-1 = 1.^- 

XO I 1 )=K0 ( i ) 

- -XOC I+0)=V0( I ) 

ia..( ( n 


c 


SET 


I N i n A L 


A I,00 ,DV V 

MAI 0o607 
MA I 0060 
_MAiaa6,c 
..MAlOUol 0 
,_MAL0Qo.l.l 
_Js!A]-LLUdI 2 
MA100613 
MA I 006 i 4 
.MA 1,0 06 1=. 
.MAl 00616 
...MA 1 0 06 1 7 
_MAi.aU,6l.o 
MA I 0 06 1 0 
MAI00620 


UP INI TIAL i’i FOk MISSION 

F-AClUi-i=;( UK + VAPa-> VP6a^‘ ( HiAPX+RLATM H ) )7 C { HAi^X+RE ART H ) =* =i« AP X-F-VP.6X ).) MAlOC 62.1 

DO .72 0 1 = 1.6 . „_MA10062.2 

,U0( I)'i=-X0ti-t-6)/Vi . . - - _MAI0Uo26 

U0t-l-+-3 )=+-X0T-i ) »i«FACTUR - MAi.0062.4 

NO=0 

COAST TUG BACK TU STakT OF FIRST uURN 


CALL-CUAST ( XQ * UO » -6UXN1 /2 • 0 • XO *uO » PH 1 • HHi ) . 
.IF.(lb.aUND.c.C*l } CU TO. .2 02 


'C-CONVCRGC-T HULC- BUR N • UU T UUUNO . -M-tS-Ci ON 
C TO INSURE CO(>lVERGuNCE 


eackwauus- in .-I 1 Mi 


MA I 0 C'026 
MAI 00626 
-,..MAL.0.062 7 
...MAI 00626 
.. MA I 0062 
-jylA iU.0600 


1 BACKS- 1 •. 

..C„. .. . 

c.SAVE Initial gtaTw uf- 

OU- 76u- . 1 = 1 ♦ 6 - - 

7S0 TUGSAVU )=X0(1 ) 

T1 MTUG=T0 

-C-PROPAGATE target - 6 1, 4 n 


ruo 


MA1006S i 
MA100632 

MAI.00666 

..MAlUOooA 
. ... MA1G0o6‘6 

MAIO,. 0.6,3c 

MA 1 Uc66 / 
MAlOOoOc 

... yAl 0063V 

. MAI 00640 
MA 1 U0o4 i 
MAl 00,6.42. 


fp FINAL bU.RN NODL - 

NL3 = -i. . 

CALL CUAST (XT . GoUM , TAUP+ 1 AUX/2. 0 . XO » UQUM , Pi t i , PH I } 

SHTT-UR— IMPULSl-VL vJ .LoCO T UN- A T. .NUDi - - .. . 

PAC TUR-= ( UK V AP Ar V f^GX-i^ ( li APX +P LAR VN ) ) / ( ( HAPX + Rc Ap T H ) ^ + 3 ^ vV AP X+ VPGX ) jMAlOOGt+U 
DO 7c0 1=1.3 MA1U0 d44 

C10( I l=XOn+cl/VAPT ; -- - - MAX0064.L 

.760. 00 ( I +3 )=~X0 I i ) 7F AC I UR MAI0U646 

C-PROPAGAlc. TU LnD of Ca--> I L'JRn MA100647 


— nq=u — 

CALL COAST (XO. GO .BURN J/2 
C SET UP TIMES Ai<!RAY 

T I MLS t 1 ) =T u 

T I MLS ( 2 ) =T 1 .MLS ( 1 ) +LjO.-:;i3 
.. Ti MLS( 3 J = r I MLS ( 2 / +COAST2 
-_T I MLS ( 4-) = T IMLS ( 6 ) +i..UisN2-.U 


O.X 0.00 .PHI .PHI ) 


T1MLS( 6) =T 1ML5 ( 4 ) +CUAST I 
Ti MES( 6 )=T i Mi:.S ( U > +r-.UANl 
C— REDUCE MASS.SLT -TlMi.N Ar^iJ TARGET ST.ATL 

T0= l IMLS( 1 J .-. - 

T.T = T I.MLST c)-- 

.. -AM 0~ ( oURN 1 +'..»U Rl -iX T I A-MDG T 


__MAX-0.Q.L.4.c 
MAI 00649 
MA 1 QOobO 
_MAi006Sl. 

... mat 0c6bL 

. MAiOOobj, 
_MAiaGCLb4 
MA I 0066b 
MAI 006bo 
.-,MAl.aQ6b?. 

MAlOOobB 
-MA1006D9 
— MA-I 0066-U 
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file: Main 


r T;i Ai'i 


CA‘.*L5i-<irjGi: .-/ionitof^ system 


MALii- -U ATL * VLLDL I -i-Y._A.NU_-UDJT-- 
AMl>DT = -A-MOOT 
VtH{ ) 

OC]W— 7^ 0 — — . ' —— 

XT ( J-3 ) =TUGSAV IJ-G) 

XT ( J )=-TULSAV C J ) -- - - - 

X-O-C^ )=..r.XO (-J ) 

> Q0(J)=-00(J) 


.MAl-w.0 ui>. 1- 
N?A I 00662 
MAI 0006 3 


202 
... 3200 


3201 


3202 


— 3203- 


XT(J-3)=TUGSAVtJ-3) ^ MAICJOOOj 

, . MA10UD66 

,.XTCJl=-rUG 6 AVtJ) M.Tm, 

-xo-c^ iAAx-u-u 

. MAIOOooo 

Q0(J)=-00(J) 

MAI 00669 

Tc -1.MPULS1 Vi.-. i - - . .. - — - , -rj 

WklTLl 1GUTPT*3200 ) .N.-UitNS 

FORMAT (//*2UX. » ^ ^•.12,*-6URN IMin’LSiVc APPROXIMATION SUMMAMA 1 00o72 

4.R.Y~-A— ,.X//-.3GX*-.*. — *-^-.5-. -.aRLii-T-AL. ..LLLMLN.f S- - * * /X *.1 2X..*. * , . 3L.Ml-MAJl-lMA.L0.0673. 

AXIS LCLLNTiv iC 1 TY , -JL 1 GH f ( APOCLL) H L I GHT ( PLK 1 GC L ) P RlOO V(APMaIOUo 74 
30GLL ) V(P£R1G£L. ) » ,/, 14X, • t 1 LOMLT LPS ) » * loX * » t KILOMETER -S) (K. I LOMM A 1 0 0 67:> 

■ ALTERS) . -.(.St: CONUS ) .. (KM/-SLC-) IXM/SLC)*! - = MAL00o7a 

WRI Tt_l iOUTPT ,320 i > A u . u. 0 ♦ HAP U , Hl>G U , T AUO , V AP O . VPG 0 ..MA10u677 

FORMAtC • T'J-^* ,r I s . 6 , r 1 3 . O . F 1 o . 3 . F 1 6 . 3 • F 1 ? . 2 . F 9 . 3 • r' 1 1 . 3 ) - 100673 

~IF-tNdURNS*LG-3>-'^Hi T.. ( 1 uOTp.T., 3 J U 2J ., .A P . LP-. HAPP.* HPGP . T A3,..* VAPP. ..V.PGP_MALL).au79 

FORMATt* PHAS ING * .F l-o.3,r 13 -3 * Fit>-3 ,r 1 6,3 *r 1 2 .2 pF9 *3 * 11-3) MAlOOot'O 

W«iTLUUUrpT,3203) A X . f:. X * HAPX * HPGX * TAUX . VAPX . ^PG X SMAlOOoSl 

—Format ( * t ransflr * *i- ip*^ *t-. 13 o*3 *f i 6 -3 ■ F 12 .2*F9,3 *- i-i.J) r-iAi ooaaL 

a'KI Tei( 10UTPT*P2U*t ) /\ T * LT » FiAP f • HPG T . T aU 1 1 VAP T • VPGT MAlOOO&J 


320.6 


.32 06 


MAI 00664 

JdA-tuu-Sau 


-MA.1.00633 

MA100669 


..3204.. FORMAT t * TARGCT * ,(-M6-P*riJ .S.Flb*3»rl6,3*l-12*2»F9.3* i ll,3t//)- . MA100 664 

WW-i TL (-1 OUT PT. *32 03 ) - Rc L 1 Hi. - - — MA-tuu-SLU 

3203 FORMAT!* RLLA’MVL ii-iCL I'NATi UN* ,Fa .3 » * ULGRErlS (McASUR-O + OR - ATMAIOOCSo 

1 PERlGt E ) * */// w:9X, *« * * # ,»S # ♦ ^ 4 # * ♦ ♦ » » MAl006ti7 

, - - : MA.1.006S3 

WR ITL ( iOUTP r , 3206 ) T • ,LG IN . o'JRN i MA1006&9 

3206 ' FORMAT!* FIRST J0;L, ulAiiTG AT*,FU-iii* SECONDS ‘ F Iv ST OURiN IS‘MA100690 

WF...1.2..2 SLCONi'JS*.) - — — - MAJ_u0.aO.l 

IF (NOURNS* CO. 2) wK 1 1 E ( I UU IPT .3207 ) COAST 1.BURN2 MAl00o9^ 

iF(NSURNS.F0.3 ) V/RlTt. ( i UU TPT . 32 Os 1 CQA ST 1 * OURN 2 , CD AST2 ,6JRN3 MAI 00693 

32.0-7— -'FUR MAT ! * CO AST - 1 6 * , F 1 0 . 2 , * SLCUNOS * , / » * SLCUNO UURN I S F3 .2.* L. SLCMA 1 0 0o9 4 

lONDS*,//) MAl00d9c 

3208 FORMAT!* 1ST COAST io*.}-i3.2** SLCur>i DS • * / , * SLCUnO OUR . . i S » * F 1 1 * 2 * M Al 0 069b 
j. SLCQNPS* */.* •- S^SONP .cUAS.T.-lS.*.*F-iu.2..-*.-SLCUNOS» ,/ . » . F iNAL-3 UR.L-U.SMja UUo9.7 

2 * ,F 12.2, * GcCONOS' ,// ) ' MAT 0 069b 

C calculate Thu WLlGurS.PA .LU ON LSTiMATLO OU«iM TIME. MAI00699 

204 iF(NUTARG-NL. I ) GO TO. 206---'-, - - MAl.OO.i'OC 

RuAD! INPU r*\A.ML.S2 ) MAI0U701 

KuL1NC = 0.0 , , ■ MA1UU702 

_2J0i> CUNTlNUu - . ■ - ;* ' — MAX.U0.7Q.3 

TOTbRN=UAl^St 1 1 MLS ( u H-T 1 i-U-S(4 ) +T l^iis ( 2 ) - MAiUO70b 

I ( TI Mils ! 6) +Ti M(-:S43 )-M r MLS ( 1 ) ) ) .MAI0070S 

__ PS 1 = thrust/! VLrR I , 1 ) -AMSOT I^T UTLRN ) - .- - - - - MAi0C706 

_ 246 1-1,6 . . .MA1GC707 

bLTA! I ) = 1 .002 MAlOOj'Ob 

lF.!a..G'I *6 ) . uLTJ-t 1 I ^-i .0 0 8...., ; — - MAL0.07C9 

248 V-eIGHT ( 1 >~3r.TA (1 ) ^^PSi/! 1 .-►ReTA.! 1 ) -IP’S I) MA100710 

203 CONTINUL MA100711 

-- ■■ ^Ai.0U7l2 

,_x .CALCUeATe LNO CUi.iO it I uNS . MAI0071..: 

CALL uVAL6(XT,U0. p-iY, TV. -1) ..MA1 00714 

..._ . ■ MAJ.(J.0.71r. 

00—449 : 1— L ,1.^ - ; — 


2Tli>- 


MAIN - 22 


FILL*: ^AlU 




i 


CAMLKIuLL MUNiTQ . SYSTLM 


241 ^ 

c * ♦ 
c 


3223 

3220 
— 32 13'- 
—32 1 0 . 
_-C TriY 


-*-C( i-)=uL>( ly — ■ 

^ '-y , 

I r ( 1 .uL *0 ) ;r ^ i 1 L ( 1 OUTP T *32 I O ) 

-IF ( 1 TUn:r4R*L0. I ) JKi I w C 1 OU fP r *32 13 > 

-IF (.1 TU.'^r*i-C.'i-:o.2 > 1 fi- t 1 OUTi’T *322U ) ■ 

-1F( ) -V4 K1 TL ( I tlU TPT * 3223 ) — 

format (* bLOlN (lU iDAVMCi.-. -Of4LY CL'h'J VLKGb'NCK • ) 

FORMAT ( ‘ AL>D 3i<0 I'URN Ai4L> RL CON VLRGL * * ) 

•FURn’iAT ( • -L'l;vjii4 f l,i ;■< i4'”^-G4i..»UNO --L Oi’w V i-R GLi4 CL « ) . . .. . 

• FORMAT ( ' f4 COPLARAk CO j4V L.-cCLNC L • • « // / ) 

TO COnVLkCL THi- CuHLkNA K . M IS 3 1 UR . IN .LL* 30 ITLRATIONS 


DO 2o0 lTck=l,30 


i-t-yi 0(2 ) +0 U ( J ) *-}i2 ) 


.247... 


-251' 


N0P=1 

,.QMAG = L>SQRT ( UO ( 1 ) 

.00 247| 1 = 1*0 
..UO ( 1 ) ='U0 I I )/UM AC 
-CALL-GUI OlK 0.0 ) ■ - 
CALL AOXUUT 
CALL CK.SFT(Ck) 

-DOOMAX=O.U 

_OTMAX-u.O ■ • 

- 00 2bl .1*1*4 

)=uo( 1 )h>uuu( 1 )-;^ck. 

T1MLlS( I )=T1MHS ( I ) +01 i MLS C 1 ) «CK 

IF < OAOS( DUO ( 1 > > * CiT • .'HjUM AX ) OOUMAX =0A BS ( 'JQO ( I > ) 

-lF.U>Af,iS( L>7 1 MLS i I > ) • CT .OTMAX) - 0 TM A X-D AS S < OT I.MLS( I J ) - - — 

,, am1=AMU , _ _ 

I F I 1 bA CR • i_u • 1 ) /(M 0 =.Vt-.'» t 1 • I ) + ( T I Ml-S ( ;:> ) —T I MLS tS)+TlMLS(4)- TIMLSt3) 

.i—4.T.i MLS(2 >“T 1 V-LS( 1 > )-iiAML>UT - - - - 

AM0-(AM0+AM1 .0 

IF ( DTM AX • LT * 1 * U-t> :|«U Ai . S ( T 1 ML S ( 6 ) ) • AND .OUOMAX.LT • .001 .AN k.'. 

• iDAbSl AMO-AMl J ,L f . V^H( 1 . 1 ) 'M.L-o ) GU TO 2S2 . 


.M^i, 0.C'7I6- 
MAI 00717 
MA 10071 rt 
.MA.LaU7.1 O 
00720 

MAI 0 0-72 I 

MAIXUL722- 

MA I 00723 
MA 100724 
. _.UAI 0 0 72 3 
MA10072O 
. . . M A I 0 0 72 / 

MA.10Q7-2b.. 

MAI 00720 
MAI 0073G 
„,_,WA10G73.1 
,.-MAl 0 0 732 
.. .MA 10.0 733 

MAI. 00.734 . 

MAI 00733 
MAI 00736 

MA1.00737.. 

. .._MAI 0,0733 

MA10U739 

MAI 0074 0. 


_.321 1 
-32 14 
— 32-21— 
2bO 
C DIO 


-3006. 
-*262 — 


FORMAT ( • CUPLaN/w^ vlic5.1U^< CQNVLROlD IN/. 13** 1 T LR AT IONS . * ) 

.. FOR MA f(* TuUN-ARUUnU, /’VCt'-ilLVLL 1N**I3** ITi-l'MlIONS**) 

-FORMAT ( * -3R3 --iiUA.'J AUUlD - AnU.. CON v'LRCLO IN*. iS.' ITLrvAri ^NS. * ) 
C O N T I i'-i U L 

not CONVFRGl: liN 3C IILKATIONS. DUMP VARlA.SLLS ANO STOP- 

- kvRlTL( lOUTPT *30U^ ) ... 

-FORMATC' PLAi'JAU ;-iiSSiUN LUO NOT Cl) Ts:.'/ LR C«- 

-Slop 

-CGNTINUL - • - - — 

IF ( I TURNk *L0 . 0 ) -Vr I T { i OU f PT *32 11) ITLK 

IF ( I TURNR . Lu. 1 ) WKl I t-: ( 1 OU TPT ,3.2 1 *+ > I TLR 

-Ir( ITUHNR.i:<J.2 ) k i ( - ( 1 OUT PT * 32 2 1 ) 1 T b R 

..-iP( I TURi-4R • bU .3 ) V.'i^ 1 I ( I GO i'PT .3222 ) • I Ti_R 


IN 30 1 TLl^AT I N 


—3222 .. format ( *- oOlUANCL- 


.y -CUNVLRGbNCc- AlHIlV, D I .N * . 1 3 . * 1 -LRATIONS. 


C ADD PLANL CHANGL 
C 

ANGLL=U . 0 

.C-.-IF- NU . PLANL CHAi-AiL SKii' AR.OUliJ 
-I F-(0 AS S ( ALL I NC ) .L T . .<b kR L‘K } -U O -1 O -20 3 


MAI 0074 I 
MAI 0074 2 
,. -MAI 00.743.. 

. . MAI 00744 
MAI 0u74b 

MA 1-0 0 7.4 6 

MAl 00747 
MAI 0074 l> 
.....MAI 0 0.74 9 
MAI U07SU 
MA1007S1 

TlAiu.07,32. 

MAI 00753 
MA I CO 754 

_MAI 0u7bb 

STOP • ) MAl 0076o 
MAI 00 737 

?4Al0O73a 

MA100759 
MAI 00760 

MA100761- 

MAl 00762 
■ ,i M AI 00 , 

.._<M A 1 0 0..7.6.4- .. 
MAl 00763 
MAl 00766 

MAl 0 0737. 

MAl 00/66 
..1A10 0 76 9 
MAiOG.7-73-. 
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I 


file: main 


FiJin ri<: Aiv 


fees' 




CA.*\b;<lL>GE MDNITQ ^ SYSTEM 


-T«AUSFO;^W-.-T*i M iSS-i O.^- l-'-SK -iNS^i. f..l ON- UF—PL A.'4L - UHAWoi 

IF I NBUKNS •N.- .o ) GLi [ Jt>40 
T6SAV=-T1MLS(0) 

T-6SAV«TlMi--S(S) - — 


-2&3- 


TlMSAV=iTT 
UQ ^ -1—1 *A 

-T-l-Md-S (~7— 1- ) - T- 1 M L 


( S-i ) 


pNAbiNG LiRbir ::nlj conditions 


Ci-vANGL IN IC DEijiiCc STEPS 


TlMES(2)-0.0 

TIMES( 1)=U-u 

N0P = 0 - • 

.-C- CALL CUIDL TO CCNLiPWL 
- CALL GUlDu(0*u) 

UO-2 5L--I— 1-.6 

254 XP HASE U )=Xt I > 

TT=Ti wLS"( o ) 

-—2540 CONTINUE 

C PERFOi^M RtliUUlRED PLANE 

RLLT= JAuS ( r<ciL i NC ) 

ANGL— DMINl t UALjS ( A.s.oLr. ) + TG • »J’TLLT ) - 

AiNGL=DSI GN ( ANGL »PEL INC) 

DANGL=.ANGL-ANGLE 

AN<>LL=ANCL 

WRI TLC lOUTPr »4C0U ) ANGLL . - 

..4000 FORMAT c ATTLNPT T J CUNVi_RGL »»F6.2i* DEGRi 

c— T-k-ST---WtIL-THLR-~-T ARGE-T---ST A rL--AT-ONL-OF-.NODLSlUNLY.- 

IF( iPHASE.NE.O ) GO TO 257 

• C RT*VT AT APOGEE OR PE R I Gf :E * RO T A T E T HROUG H ANGLE 

L)0--?5G.- 1 = 1 .5 - — -; — 

256 XTt I )=RT ( i ) 

XT( 4 )=C.O ■ 

{ ) =VTC2 ) aOCUG ( AncLu*OLGCON) , 

XT (6 ) = \/T ( 2 )4DS IN I ANoLL^J^Of.GCQN) 

GQ TO iLiOo 

c-check- Tu--sl:l; if. -convi'.kg ing - uacxw ards - 

...257 IF ( luACK .LC . 0 ) CD TO 2o7o 
. -,-C. CQNVLnG i NG BACKWARDS RUl'ATL XC AND DO 

CALL- - ROTA T li ( XO *Priot'ir * OA i4k-*L4iO;_GCuN.* o ) — • • 

CALL ROTATE ( OC *PRGi,.L , DANGL-^OEGCON * 0 ) 

GO TO 25v 

C- COtHVERGi nG - r DR .V AN JS * F 1 NtJ Ir . 2 .. OR-. J bOMNS ■ 

__2570.-- IF ( NbURNS • L.^- «S ) GO I Lj 2 Stf 
. _.C.-.2~BUKt4 M iSS 1 urv #ROT AT rs<i..SDi-H xT 

CAU RUT AT.LCxT-*PKG{-U,UA,vGt-?i'DLGCUN ,6 ) . . ... 


AiO.U 7.7.1. - 
MAI 00772 
MAI 00773 
- MA. 1 0 0.7.7. 4_.- 
MAI 0o77b 
MAI 0077o 
-_MAl.auX77- 
MAl 00776 
MA I 0 0779 
.-,/.A.l-0 u76a - 
MAI007S1 
MA100762 
-MAI Q0ZS3- 


MAI 00784 
MAI 00785 
. MAI-00736 
. MAI 00787 
..MAI U0V68 

_MAxoo.7a.y.- 

MAIC0790 

MAI00791 


GU TU loco 

(2 3— burn forward i^UTATL XPPJASr. 

^258 .CALL -R l.l'l AT rl ( XPHA St- ♦ t ’R GL L, * DANGl ♦DL v»CiJ N* b ) 

--C...SET UP 3-3URN AS 2-GURi'j 
,259 DU 2o00 i~l*D 

-.2600— -XT.l-1 )-XPHASL I i-> - 

C ON FIRST f^ASS Sh.T Ui' AN INITIAL COAST 
1F( ISECD.LO* I ) GO TO lOgO 
ISECD=l . — - - ■ 


■ NU-0 

TbCT=-~50C. 

C ALI. -tGRA-S'l ( Xu » wG *-Ti i .* Xu *.U0 »-Pf t I ♦ Pf-1 1-)' 


'• “ " ' * 

' 

, MA1U0V93 

PLANE 

CM NGE • ) 

MAI 00,794 

PMA5L 

...NOg. CALLLO-.)- 

-_i.lAJL,0U7Ji5-. 
MA 10079b 
MAI 00797 
MAI 0-0-7-95. 



MAl00/'99 
MA XOObOO 
_MAl.0.Ub.0,X.. 



MAI 00802 
MAI O0S03 
..MAl-UU5u4 ' 



MAI, 00805 
MAI OoaOb 
MAi.uoa.o7:.. 



MAI 00608 
MAI 00809 
.J^Al 008 Xu 



.MAI 0081 1 
MAI 0 ua 1 2 
MAI 0 06 13 . 



MA1C0814 
MAI 0081 5 
.MAX008X6 



M A I 0 0 a 1 7 
MAI 0081 a 
, MA-IJJ..O£ii 9 



MAI 00820 
MAI 0082 1 
MAi.O0.82Z.. 



. MAI 0 082 3 



U-1A I 00824 



MA-X,.G.U8Z5- 
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file; MAIf^ 


TH A,-«l f-*l 


CAMt3RiDC,Li MONiTOr- SYSTEM 


0 1 {.) C-T - — . - - • 

C“ REDUCE INITIAL MASS TO .<r:rLf.Ct 
IF ( ibACK .Lu . 0 ) Oil tii luOO 

AMSAV-AMDOT*t T oSA V -T ;-.SA V ) 

VEH( 1.* 1 ).=VLH( 1 * 1 )+AMSAV . — 


ESTlMArED third bURN 


..MA1U0L26. 
MAlOCjii^y 
MAI GG62ii 
-^5AL0 06^'j.. 


lOOU CQwTlfiUL - 

Gr-SE-T~UP— i_ND - CUND I T l UNS -.r.b;<-.'-EANL 

call ti VAL6 CXT,bO * PTv/ * tv , “1 ) 

DO 261 1=1 *6 

—26^1—- -CC( 1 ) =OD ( 1 ) ■■ 

C - - 


30 ITERATIONS. 


C TRY TO CONVERGii WITH PlAjmL CHANGE IN .LE 
DO- 29.0 -_i.TER = l • oO ... . . 

UMAG = DSQR r (QO I 1 ) i E)=4==J‘2 -f-OO ( 3 )'T4-2>. . . 

DU .209 i = l*u 

.2-^9 Ou^-l3-OU( L )/uMAG ■ 

CALL CUl Du I 0.0 ) 

CALL AUXUUT 

CALL- CK.SUT ( CK ) - . 

UQU.MAX = 0 . V 

UTMAX=0.0 i 

j„PHi.NT-=.0-- 

DO 2^3 1=1.6 

DO ( I )=00 ( 1 ) + D00{ 1 ) =4--Li< 

TlMLSt 1 ) =T IMES{ U + L.T 1 MuS( 1 )*C.< - 

IF C DAOS ( UQU t 1 ) ) . C I . -JG-.A/. ) --DU 0 ( 1 ) 

.283 . I F I DAbS U->T I MES ( i ) ) . o i .UTmAA) DT maX= 0 AbS ( L'T I MuS ( 1 ) ) 


MAIOObba- 

MAl-0 0631- 

M/LlOLAiSE^- 

MAI00633 
MAI00634 
.MAI.OOtSSS. . 
.MA 1 0 0836 
MAI 0083 7 
-MAl.aQ,a3B— 
MA10O839 
MAI 00840 
..MA.LC084.1..- 
.MAi00S42. 
Mz\l 0 064 3 
_MA10.Qb.4.4„ 


1F(13AC<*L0.1) AM0=-Vu; « ( I « 1 ) + ( TI ML S (6 ) -T A MLS (6> + TIMES{4) 
I +T iMLS( 2 >-”T iMLS ( 1 ) > -AMuuT 

AM0 = ( AMl+A.’lO J .6 

I F { D TM AX • L f • 1 • D—o ^'‘lA..;S i } 1 ’’iLS ( 6 > ) .AND .DOOM AX •LT..00I.At-f 

l.,..UAu S(AMi“.AMO)*LT.Vti»'iCi«lJ^l.D— b) oU T(J 1 


ITLRAriONS. DUMP VARIABLES AND STOP 


-29-0— CUNT I NUL 

C DID NOT converge IN 30 
WRITER 1 OUT PT .3 06 6 ) 

.3006 format (* cjut-uf-pla.n^ m I ss ion - d id . »\<jt converge. In 30.-I 

ITGP . * ) 

STOP . • 

WKl TE( 10UTPT.3212 J A Nt:-L L . I TF. R 

3212 FORMA T ( F7 • c * * DEGRiL PLANi_ CHANGE CONVERGED IN*»I3.* 1 

IF ( DAOS ( Ai-JGLL) .LT.i.> AwS( RELINC). ) - GU TO 263 . 

.C-- CHANGE LACK TO J— 3URU ir NLCESSARY. .. . 

IF ( NbURNS . NE.6 1 OJ Tp 266 ..... 


MAI 60845 
MAI00846 
.MAl.a.0,b,4.7_,. 

MA100O4E 

MAl0e64D 

MALQ auO.0.. 

.. MAIUuabl 
MAI 0 0652 

... MAIUU653.. 

MA106634 

MA100655 

dAl,0.,0b56., 

-TIMES (3) MAlOOSbT 
MAI 00656 

MA 10.06 59 

L>. , MAI 00660 

^ .MAIOOObl 

MAl.O 0.662_. 

MA100H63 

MAI 0 0664 . 

THRATl DNS. _SMA1.6 06.65. 

MA100666 
,. MA 100867 
M AIOGSEl:^ 


MA100669 
DERATIONS. * )MA1 00670 

MA 10087,1 

. MAl 00872. 

. MAI 008.73 


292 

— .-UU — -i -A— ‘+— - - - — — 

T 1 MES( I ) =T IMES ( 1 +2 ) 

T I MES ( 5 > -=TdSA V 




MAI 008 75 
MAl 00676 
MAl 0-0 8.7 7 






MAI 00876 

, 

XT ( 1 +3 ) =“TUGS/\V C i +3 ) - - 




J« 1 A 10 08 79 

— 293 — 

— XT (-1 > = TUGSAV( i > r — '4 — 

........ 

: 


MA -1 E 0 . 6 EO.. 


MAiN 


FILC; MAIU 


p iJi5Ti-< A IS 


I 


CA^■.3R ibcVil MONITOR SYSTEM 


^ T-T;=.TiMSAV 

" «ELlNC=O.U 

C ADD mass to irvllTiAL ma:- 

I . I-p.( -I Li A CK-* tiS* 

1TURNR=2 

X-ROTATE -X.T....IF. FORWARDS MIS 


IF liACKWAROS MISSION 
.1 ) = l ) -AMSAV. 


TON 


I-FT-i-U A CX O -i )— DO— T-U-.-2 DO- 
DO 294 I = t * S 
X7 < I )^RTC I ) 

— 294- XT( I-1-3 )-VT ( 1 ) -- 

CALL k'OTATcCXT .i‘UOLL. A 
CO TO 2C3 

C--T &-S-T— Tu SOL . IH. -:*USSUJN .. rURNAUOUND - 1-1 
295 CONTINUE 

IF( UiACK-NL. I ) CO TO 2o0 


.TslAJjTv) EE-1.— 
MAicoaoa 
MAI ODiiOS 
-JAAT O 0aS4. - 

MAlOOoSS 
...MAiooaao 
M A I o ci.Ea -Z — 



MAI ooasa 

MA I 0Dti39 

u^A-i.axia90- 

. MAI UUE9 1 
MAI0U692 
^ilAXaGE-y-i.. 


ilOiM Aris^ i\ _,CUNVLriGE. 


,.C TURNAROUND MIS 

-,C 

_C—SL-T— UP— times- ARRAY-»00 .AMD -AND xC 

WRITE! lUUTPr.AoOl 1 

4001 FORMAT! • ruRNAKOuND MISSION*) 

-P t s J-I Mils ( S ) “T i i‘-LiS ( D ) — - 

C1 = T IMeS( S 1 -TI MLStA > 

b2=TlMLS(4 )-l 1 HuS ( s ) 

C2-=I-1MLS(S ) I MLS ID ) - 

03=sT IMLS(2 ) -T I Mi_S ( J > 
TAL1GN-TT-T1ME.S( 6 ) 

-T TLONsT I Mi.S ( u ) — TO - 

T0» fALlGN+llMTUG 
TIMES! 1 ) -TU 

TlMtiS! 2.)=T1 MES < 1 i 


T 1 MLS {3)=T'rMES‘(Z)'+"CI ' 

T1McS( 4)=^TIM,-:S!.:' }+iC^ 

— TIMES! 6) =T IMCS ( '+ 1 +CL 
TTMLSt b ) =T I MLS C ) +i.o 
tt-to+ttlgn 

_A.MOur = -AMDOi' - • 

VLH < 1 * 2 ) ! i . L ) 

DO jlO 1=1*S 

-Q0(l)=Utl) 

QO ( 1F3 )--0 ! 1+3 ) 

AT! 1 ) =X0( 1 J 

..-XTT.1+3 )=-Xu I 1+0) 

aMQ=VEH( 1,1) 

C COAST TUOSAV TO ALioN v; I j M T 1 ;l%S ( I ) 

N0--1 - ' 

C A L D C OA S T ( T UC S .A V » 3 
C. KLCONVERoE with r-UKWAi<L.> 

IiiA-CtC-0 . — 

Rl.L1NC=0«0 

ITURNR=1 

GO-TU- 203 


-31X- 


UM , T A L IG N * AO • ODUM • PH I * Phi 1 ) 
•’IS a I UN 


MAI 00694 
MAIOOB95 
„_JWA1.0,Qb9.S- 
.. MAi 0039/ 

MAI 00393 

MA I !; U39 .9., 


...260 CUNT INUE 

._C--VE.RiFY final ORbIT. 
DO-2 o2- -I =1 .3 


MAI00900 
MA10090 1 

MAi-0.U.902,.. 

. MA10090D 
... .MA100904 

MALU 09 .U. 5 ... 

MAI00906 
MAI 00907 

MALOaS.GS 

MAi 00909 
MAluOVlO 

MALO. 0 . 9.11 

MAI-D0 91.2 

MAI00913 

MAL00914 

MA 1 0u9 I D 
MAI ooyio 

MA10.0 91J'. 

MAI 0091 a 
MAi 00919 
MAI 0 0,92 0 
MA100921 
MAI 00922 
_MA.i,0 0923,- 
MAi 00924 
MAI 00923 
..MA100923 
MAI 0092 7 
MAI 00923 
_MAl-0.0,y.29-. 
MAI 00930 
MAI00931 
.J^A 1.00932 
. MAI 00933 
JviAI 00934 
Q.93o, 


F 1 Lc I A I iN 


I • i.j K r i \ ' 


I 


CAv.iJr'^ .MUNiTO 


.YSTLM 


262 


— RTA( i )-X( i ) 

VTA ( I ) =X( i +3 ) 

Call Li_MjNjT s ( r< r a * v i a » l a * h a * pc a » t au a ) 


WPITL (-10LTPT.*3^i 3 ) i A * V f A • A A . LA * HA ,PGA *TAUA 

■ 32 Fl)pMAT ( 1 UKLiT- ACTUALLY. ACHltZVLO; ' - - 

.I.« PQSi TI JN - **301 A. O*/.* VLLOCl 1Y-».301A.6*/.- 

-2-'LuiLM-I--MAUi3p--AX-I.S=-*.»r -La*2., / .-LCClKT.P I Ci'T Y = ? G* 6. 


MAI 00937 
MAI 0093L 
_MA1.0 0.93 9 


/ 


3 • H-VLCT0P=* •3l>i “» -o, / ♦ • 

BUP 1= r IMLS ( 2 )-“1 I -=^S( I > 

— = T I MLS ( a ) — 'i i -YL;} ( o. > 

. .bUP 3= riMLS(6)-TI 'i cS ( ^ ) 

.... C0AS0=TiML3( 3) -T c 
— CUAS4~Tl.Mr.SM ) -TO- -- 


PuP 1 GLL= • *3 l)14 .6*/ • * PLR 10D=* * . 1 0 ,2) 


MAI.00940 

MA1U094 1 - 

MAJLaoa.43- 

MAI 00943 
MAI 0094a 

.MAl.Q 09.43- 

_ MM. 00940 

MAI0094/ 

MAi-OO 9 . 45 - 


CQAS2*TIMES<3) -T1ML3(2 ) 

C0AS3=TlMLSlD)-f lYh0t4) 

XF(.NoUPNS.lU, 2 ) ..vVPl fL( UJUTPT *32 1 7 ) CO ASO . 3UK2 t CUAS3 •OUR 

32J.-7..-. F03MA n//» ' CGOVLi^LLu CUASTL A^D LUPTmO 2-liUKN MISS 

1 initial CQAS i ~ * *r lo - 2 . F I ST DUPN = ■ i F 1 0 . 2 • / * * 

2-=-.»-*F-iO*-2* /*-* .F i-NAu-.-L.Ui<N.-l * F-1 u .2 ) 

C 

) wPI Ti.; ( 1 OUTPl *32 13 ) CQASi *3UR 1 * CU AS2 *BU 


IF ( NOVii^NS »LQ • 

32_lb..— .FUWMAT t * UJUVLkcuU CDaSTS AND bUUr4S FUP 3-3UPN MISSION 

1*., initial CUAST= * .r I 0 ..3 ./* * FInST 0UKN=* iFI0*2i/* * 

2 * I F 1 *2 */-*•••- SLCLNP jjUPN~*-.*F- 1 0 *2 * - CUAST.=..* *F 

3» FINAL HUKN- • » r 10 . 2 ) 

C 

— C-GUIDAnCE:- SLCT ION 

...c.. 

C- . 

— C-CHLCK- I F--GO-1U ANCL -DU.S-li'H.u-*- 


MAi 00949 
MAI 00950 

MAI. 0 U 9 O.I 

.!UN: ■,./» . MA10095Z 

; LCGND.-.CQASTMAI 009S3 

aAi-OU 9.3.4... 

WAl 00955 
tCOASS * bURM AI 00935 

WAI.00,957. 

;• ,/* MAI0U955 


IF INUTARG.cO*-! > STOP 
C CHECK IF 2 UR 3 BURNS. 

IF (.NbUf^NS.i.:0.2 ) CO (H 4X0.- - ... 

_,C-3.-bURN MiSSiUN. RLf''uVi.. ANY I N 1 T 1 AL COA-ST. 

IF ( OAuSI T U-TTMLS ( 1 i ; . LT . TlRRUR) CU TO 4o5 


CALL COAST ( XO * 00 ♦ r i :ylS( 1 )-TU .XO *0U.*PHI •PHi ) 
T0=TiMLS( 1 ) 

.406 AM0 = VEH( 1*1) - 

^ TOINT- TlMiiSI 1 ) 

IF ( 130UNO .EC . X ) on Tti 4 10 

4>-,TURN -AR0UNL> OU rbOUNU .. 3 -oOkN ...M-lSS- iUT*., 

c ■ , , 

C SET UP STATE AND TlMu AP.PAYS. 

L>(j i =. 1 , 3 

XTS I i ) =X0 ( I } 

XTS ( 1+3 )=X0 ( i + J ) 

OTS(-I ) =00 ( X ) 

OTS( l+3)=-O0 ( 1+3 ) 

AOC 1+3 )=-X( 1+3 ) 

— UO I I +3 ) =~0 C 1+3 ) 

XT (-1+3 )=-X f3( I +o ) 

— : _..xo ( 1 ) = x c I ) 

-oo{i-)=c(i) - - 


• tLCOND COAS.T.MAIOU959 

MAI JjJ 9 du.. 

MAi 00961 
MAI 00962 
-.MAI 00i)c»3 
MAI 00904 
.MA100905 

M A lX).C. 9 oo. 

MAI 0096 / 
MAX 00965 
.— MAX 00.969 
. MAI 00970 
MAI 0097 1 
_MALc u972- 
MAl 00973 
MAX 00974 

MAi 00975 

. MAlU097o 
.MAX 00977 
— KALa097L . 
MAI 00979 
MAX 00980 

MA100981. 

. MAi0o9bL 
. MA100963 
.-jAAlUi3964.. 


MA I 00955 
MAl0w986 

MAlCLQua./.. 

. . MAI 00953 
. . . . A 1 0 0 9 5 9 
_MAl0099o 



MAIN - li 


file; main 


, iM 


CAMbr^lDGC MONITO.^ SYSTEM 


406- 


_.X r (-i-)-XTt>( I ) 

til=TlMLS(£)-Ti.V:(_S( I ) 

Cl=TlMtS( J ) -n MI_S (2 ) 

-as^T jmls( A ) -ti-ml.v - 

C2=T 1 MES ( b ) -T I t ^ ) 

. bi=TlMuSt6)-riMi,S(;3) 

— T-U4t;S-(--l-)-200G. 

TlMtS( 2 > = T IMES ' 1 ) 
times( 3) =T1V^lS ( •? i +c;!f 

-TlMEStA>=TlMESt,>)+..A - — - - 

TlMES(b)= riMES(4 )+C 1 
T IMES( b) — r i^cb ( 1 

TT=T1MES(6) 

BURN r=DABS ( TIMES ( 6 ) - I 1 MEb I b) +T I MESl 4. )— TIMES ( 3 > + 

l-TlMtS(-iJ)-TiMES( 1 ) ) — - 

AMO=VEH{ 1 . 1 ) -Li'JNo^T'J'VEri ( J ,2 ) 

AMDOr~-AMUUT 

VEHt^l ..2 )=-VEH( 1 ,2 ) • — , 

SET UP INITIAL COAST FOR uACKWARDS MISSION. 

NUbO 

CALL- COAST (XG*00« -bO k, • . X u *CiU-.P»n 1 *PHl ) 

_ . TO = TU-buO. 

.SET UP MONTL CAkLU NUNS. 


_MAi_UOaO-l- 

MA100992 

MA100993 

MAT.LU994-- 

MAI0u99b 
--MAI00996 
—MAI 0-099 7— 
MA I 00996 
MAI 00999 
.—MAI.Q lOGO. 
MAIOlOOl 
MA101002 
— MAIU1O03 . 
MAIO 1004 
MAIOlOOb 
_ MAI-01006 
. M A 1 o I U U 7 
—MAIOIOOL 
MAIO-1 0U9 . 


. C 

c 

- C CALCULATE LNi> CnNOITlONS 

— CALL" OVALS I XT ♦ OOUM ♦ P r V • r V* -1 ) — - 

UU 4ll l=l»o 

.411. ecu )-U0( X > 

— G-a A V S-:l |4 1 -T-i . A L; , c 0 IS* 0 i TlON L u.< i.LXT..Mu.4TL- CARLO - A UN 
DO 41b 1-1*0 
XoS< 1 )=XGl I ) 

oosri )»wo( 1 ) - 

V£HS( I )=V|-.H( 1.1) . ■ 

TiMLSStI J-TlMLSt 1 ) 

XTS-( i ) =XT( I ) 

cese I )=CH i J 

VLHS(7 1.7) 

-TOS=T 0 - - - . 

TTS = ir 

AMoS=AMO 

— 1 hOUNS-1 uOUNO - — 

C LOOP FOR MONTL CARLO iXUiMS. 

IPRINT=1 

vDO- 4?0 MUNT*_-=i .MCA;sLO •--- 

__ — . XUPOAT=G 

e.RESTOR^^ VANlAbLuS. 

• - - oQ .._4iib • i =• 1 . t> 


415 


VEH< 1 . 1 )=VLHS( I ) 
xo( I > = XOSt 1 ) 

~U0-( I )-OGS (X ) 

XT( 1 > = XTSI 1 ) 


. QT I 1 > =OTS ( X ) 

r:i.ML-b( li^TlMLG 


T i J. 


MAlOXOlo 

MAI01017 

-MALUlUXS 

MA101019 
MAI 01020 
-_MAl.alXJLl .. 
“MA“I‘0'ia22 
MAI0102J 
..MAl U L024 
MAi0l02b 
MAI01026 
_MAl.0.i.02.7... 
MAI01026 
MA101029 
_MA10L03g 
Max 0 1 031 
MA101G32 
_MAi.0-l-U33 
MAIO 1034 
MAI0103S 
„X.1AI0 L036 
MA101037 
MAIO I 036 
_MAi-ai-03L9 
MAIO1040 
MAIO 104 1 
_MAX01w42 
MAI01043 
. -MAXO 1 044 
MAXL 1 04-^ 


MAI 010X0 
MAIOlOl I 

t^Al,0.LO12 . 

MA101013 
. .MAl 01014 
,'4AUXlCUb- 




MAIN FUiiT:^AN i-M 

CCt-i )-CCSt I ) ■ - 

Vcim i • 7>=^VnHS( /) 

DD 41^6 I =c; * 1 

-DQ-4-2o-J=l *7 ' - 

-ViiH{ I * J)~O.U - - 

-lBQUNL)-=lv>UL»i-ii> ■—; 

»-i ) 

TCUOCK-O.O 

TACCUM-O.O 

-AMO=AMOS- ••• •'• — 

..TO = TOS 
-TT=TTS , 

-lF-( NiJURNS *3 ) CALL— iiCwCii(4 bOUi-JL/* TO 
IFCNRbRNS.FO.S ) CALL CuCb 
CALL STATir>(MCARLO> 

-STOP 

• cNU .. - . 
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INT } 


MAll 

MAI ( 
MAI ( 

MAi I 

-...MAH 
_.MAI I 

MAl-i 

MAI ' 
MAH 
. -MAli 

. ..MAI' 

MAi- 

MAl' 
MAI 
_ , JrtAl 
MAI 








AUXOUT 


Subroutine AUXOUT 

A. Purpose 

AUXOUT prints the status of the convergence, from the most recent 
can to GUIDE. 

B. Input/Output Definition 


Input Parameter 

Symbol 

Definition ^ 

VD 

II 



X 

X 

Vehicle final state 

XTF(I) 1=1,6 

^T 

Target state at same time as above 

TIMES(I) 1=1,6 


Array of times at ends of coast and 
burn arcs 

Q0(I) 1=1,6 

^0 

Costate at start of mission 

DTIMES(I) 1=1,6 

hi 

Requested corrections to TIMES 

Dq0(i) 1=1,6 


Requested corrections to costate Q0 

lOUTPT 

- 

Output device number 

Output Parameter 




None. 

C. Method of Computation 

The only variable calculated is the estimate of the total burn remaining 

COST = |(TIMES(2) - TIMES(1))+ (TIMES{4) - TIMES(3)) 

+ (TIMES(6) - TIMES(5))| 
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CAMijf^lDGL MU.\JTO): SYSTcM 


— SUtii?OUTi.Ni»_--AUXOU r 

implicit RisAL’.'C^C ) 

CDMMOM /GI D I N/ A T Co ) . i T * X -j ( o ) . TO . A MO , Vc H ( 10*7) ,Q0 (O ) ,T 1 

— LOMMOK -/OOLVIC /IOUTP f 

common /o I OOU r/L)viOCo;*oflMi-G(o),u(l<i*12)*DC(12)»X<6)*Y 
1ZU2* i 2 ) ,0 ( o) . OOMMC ^ ) .OM 

WRITEC lOUTPT*! IX.XTF 

FORMAT ( // * * X c O' »TA I Nt:L> > = • *6Li iA.o*/** X(i;i£SIRLD)='*6fl 
-COS T =0 AuS ( -T I G t o ) -T l.ML.S(S) + TIMi:S(4)“T IMtS ( 3 ) + Ti'MLS ( 2-) 
WRITE! lOOTP T*2 ) oUliT 
FCjRMATC* UEMAiNli>!o _,URNi- • . O 1 4 . o ) 

— WR l-TL-t-l OOTP-T-*P )- O U ..OO O MES '• 

FORMAT (IX** 00 - * . C.i: 1 < > . o * / * * DO 0 = • . 6i* 1 6 • « * / * 6X » * T = * * 
lbL16»t>*/**+X»*Dr— 


A.UX-0 0 UU-i.- 

AUXOOC’02 
LS< 6) *CtO> AUXOOOOO 


AUX00006 


AUXOOUOB 
b) AUX00009 

-T-IM LS U-). ) AUXU.UL I O- 

AUXOUU i 1 
AUXUOO 12 

Auxo 0 ai .0 

AUX 00 0 14^ 
AuxbboitT 

^.AUX ouo.l 6 

AUXQ0C17 
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Subroutine BCBCB 
A. Purpose 

Subroutine BCBCB is used during guidance mode to take the vehicle 
through the first burn of a 3-burn mission. It operates in either a back- 
wards mode (outbound mission) or a normal moc^e, and is called by MAIN at 
the start of each Monte-Carlo run. It in turn calls FORWRD at regular In- 
tervals until the end of the first burn, at which time it changes mode (if 
backwards) to the normal mode and calls CBCB to handle the remaining coasts 
and burns. BCBCB also modifies the TIMES array on each cycle to reflect the 
fact that part of the first burn has occurred, calls GUIDE to reconverge the 
mission with the new (possibly perturbed) vehicle state, and adds the resultihg 
corrections to the TIMES array and costate. On the indicated cycles (lOUT = 1 


or next-to-last cycle in the 

burn arc), 

subroutine NAVOUT is called to collect 

the Monte-Carlo statistics. 

On the last cycle in the burn arc, the call to 

GUIDE (and the addition of the corrections to TIMES and Q0) is skipped and 

CBCB is called with an initial step time 

of zero. 

B, Input/Output Definition 

Input Parameter 

Symbol 

Definition 

IBOUND 

- ' 

0 - outbound mission (implies backwards 

mode) 

1 - inbound mission 

T0INT 

- 

In backwards mode, the actual value 
of T0 

TRUEMS 

- 

Vehicle mass before start of burn 
(normally equivalent to AM0 except 
when in backwards mode) 

XT(I) 1=1,6 

^T 

Vehicle state in backwards mode 

TT 


Tine at start of first burn in ojtbound 
case 

T0 

^0 

Time at start of first burn in inbound 
case 

lOUTPT 

— 

Output device number 
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Symbol Definftion 

.1=1; normal guidance step size during 
burn 

1=2; not used in BCBCB 

- Vector of times at end of each leg (or 

start of each leg in backwards mode) 

Output Parameter Symbol Definition 

IPRINT 

MODE 

AH0 

TIMES(I) 1=1,6 

Q!3(I) 1=1,6 qQ 

T0 
TT 

CC(I) 1=1,6 

C. Method of Computation 

After zeroing the time accumulator (used to determine when Monte-Carlo 
statistics are to be collected), saving the vehicle initial mass, and ini- 
tializing several control integers., BCBCB branches to one of two separate sec- 
tions of code, depending on whether a normal 3-burn mission is being run. In 
either case, it is assumed that the first burn begins immediately, with no 
initial coast. . 


Always 0; shuts off printout resulting 
from calls to GUIDE after first step in 
first Monte-Carlo run 

Always 0; restores mission to free-time 
rendezvous (backwards mode only) 


Mass at end of first burn 

Vector of times at end of arcs, with 
first burn deleted from the vector 
(TIMES(1)=0, TIMES(2)=0) and, in back- 
wards mode, the- vector restored- to its 
normal form 

Costate at end of first burn 

Time at end of first burn 

Time for which target state is valid 

New end conditions for target (backwards 
mode only) 


Input Parameter 
DTYPE(I) 1=1,2 

TIMES(I) 1=1,6 
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In the backwards mode, the TIMES array as supplied to BCBCB is already 
reversed and ready to use, as are T0 and TT. The weights are set to 1.0 since 
the backwards mode works best with hard constraints and mode is set to 3 to 
change to a fixed time rendezvous. Subroutine FORWRD is then called every 
DTYPE(l) seconds during the first burn, with the exception of the last two 
steps which are approximately equal to ei.ph other and less than DTYPEt'l )./2, 
and TIMES(6) is updated. Each time the print accumulator exceeds PTB, subrou- 
tine NAVOUT is called to collect Monte Carlo statistics, and the accumulator 
is reset to zero. M/C statistics are also collected on the next-to-last step 
in the burn arc. Also, following each call to FORWRD, except the last, GUIDE 
is called and the corrections are added to Q0 and TIMES, Q0 is maintained at 
unit magnitude, and the estimate of vehicle final mass is recalculated from the 
mass rate, current mass, and requested changes in the burn times. On the next- 
to-last call to FORWRD (NLAST=1), subroutine GUIDE is called repeatedly (with 
no changes in vehicle state) until the miss in final position is less than 1 
kilometer. On the last call to FORWRD, GUIDE is called but no changes are per- 
mitted in the TIMES array and Q0 and the weights are restored to their original 
value. In addition, the flag is set to add the navigation update corrections 
to vehicle state on the very first call to FORWRD from CBCB. The mission is 
_then„turned. around to normal mode, and t h e ta rget end co ndit ions r eevaluated. 
Finally, subroutine CBCB is called to handle the remaining coasts and burns. 

In normal mode, BCBCB works in much the same way, except that the states 
and TIMES array are not reversed, and T0 is updated rather than TIMES(6). 


BCBCB > 4 


ZERO ACCUMULATOR (FOR PRINTOUT) 
STORE INITIAL MASS 
SET NOMINAL STEP SIZE (DTYPE(l)) 
STORE TIME AT END OF FIRST BURN 
(BACKWARDS MODE ONLY) 

SET STEP INDICATOR VARIABLES 


1 


SAVE WEIGHTS; I 
FORCE HARD ic — 

CONSTRAINTS 
(WEIGHTS=1.0) 


YES 





X BACKWARDS 
< MODE 

(IB0UND=0^--^ 


SET STEP 
SIZE FOR 
THIS STEP 

_U 

CALL FORWRD,' 
INCREMENT 
I TIMES(6) 

L 


-LAS.r^ 

STEP IN \ yes 
BURN 
•EG. 



CALL GUIDE 
AND CKSET, 

ADD CORRECTIONS 
TO TIMES AND Q0 


CALL GUIDE 
TURN MISSION 

AROUND;. -SET TO 

SOFT CONSTRAINTS 
(WEIGHTS?<1.0); 
RESET INITIAL MASS 
AM0; REEVALUATE 
END CONDITIONS 
SET NAVIGATIONAL 
UPDATE FLAG 




COLLECT 

M-C 

STATISTICS 




CALL CBCB 
AND RETURN 


i 


CALL GUIDE. 

CALL CKSET 
ADD CORRECTIONS 
TO TIMES AND Q0 


N=N+1 

MISS - END 

CONDITION 

ERROR 












n o 
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FILh: bC^iCfi 


r»\ AU 1 ' I 


CA*-<BKID(jL MONITO SYSTtiM 


-bUUV^UUl-IWti- 


aCtiOiV-O-UJ. 


.jiCitC.;; ( ♦-tU-INT-) — 

C SUHKOUTINE£ m TAKt fHt: TUCj THI-iOUGH TrlL INITIAL HURN OF A 3-QURN 

C MISSION IN GUIOANCL WOOL* /;0;<NS F05< BOTH 1 N»UUNO ANlj UUTliOUNL) MI SS I ONS BCB 00 G 

X«PL IC-J T-RLAL*^- i A-H»-0-.-i ) 

COMMON /BVLUU T/Ar ( o ) 

„ Common /uoDATt^/iuPoAr — 

GOWMOti— /CPHY-S/-Ulv.» :<t„AN-T-j-i-*UOM-l-*;)UM2..»£iUMO — 


BCb0002o 


4;lCaQ-Q.a4 0 
bCBOOOSO 
OCOOOObO 
UCaOXLi-7.U- 


SCBOOOFO 

BCB00090 

_43caaai-a.o 


COMMON /ONL I ! lE / i Pk i ( J T 
COMMON /CMU'Ji./MUOd 

■COMMON- /CNAV/1 rOLMS-* i'NA V , TACCOM » I got 

/G i o i rj/x r (, o ) » r T . xc i u ) , ru » amo * vlh ( 10 * 7 ) .uli ( 6 J *t i mls ( c> *cc ( o ) bcooo 1 1 o 
/ u XOLUT/ijO u ( 6 ) , u'l i MLS ( 0 ) -t L ( 12 * I O ) *LiC ( I2J «X(o) *0-b)'pZ( i 2 » 12 ) 0 CB 001^-0 
SM£o )- 


COMMON 

Common 

) (-S ) 


.kiLrlO-CLL-i-O 


COMM-UN. -/C INOLX /MARC * I. Ai'i C * JMA X * J M * JM AX 1 * BLAST » NO * NOP »NF< ! GOS 
COMMON /CPCNT/PLRT ( 3 ) ,uTrPLl{ 2 ) ,oFAC 

COMMON- /CJ/aFTA(o) .P;,I *UF • CK. * Ti uJUnO * UaOUNt> 

COMMON /C V/ r/v.’tl ia^^T ( i;> > 

COMihON /C\JCST/oI Co) 

-- Common - /-ao a vi;-/ x i s c t> ) ♦ o-t s i-g ) — - • — •— •— ■“ 

DIMCNSION WTSCO) 

C ZERO TIMC accumulator. 

TACCUMsO.O- 

.. , lOU fab • • 

._..C..SAVU INITIAL MASS 
TUT LM& 3 .TROLMS 


BCBOO 140 
bCBOO ISO 
BCtiOalbO- 


- LiCBOO 170 

uCBOO 180 

b-CB aO- 1-9 0 


BC000200 
BCii002 1 0 
B.CaOU2ZO- 


,:.iJCb00230 
. BCB0024U 

.^ct>aa2-BG 


c SET nominal step SIZE IN bURN 

DT=0TYHC( 1 ) 

-C—SA VL-T 1 Ml • at C iT(.) OF 
IF I IbQUNO.tiu.U > 

C SET VARiAoLLS 

NLAST-^ 0 - 

last=o 

NPOINT=0 

lF-( 1 FOUND. C:CA*1 J 


last •i..UWN -(=TlMri AT- START - UF-.-1 ST BURN- • ^ ORW AN DS>L 
T u L .4 V L- T .1 Ml S t O ) 


BC300260 

BCB00270 

^CLiO.02.80. 


aCB00290 

bCbOOZOO 

^CbaubLO. 


GO T. 


i oo - 


OUT liOiJi'4D- M I S'_> i 3N * • NO v;- . . i. F4 --GaC.^ W AR O.'-v - • •• - 

c 

c 

C--St£ T - Vi/F 1 GH r a . • * 

C . WUGHT { 1 ) - 1 RLFlLCTC UAi<L> C O i'.'SlK A i N T S 014 BACKWARDS BUlTN 

DU 1 1 -1 * J 

1 ) =WF I GH J ( i ) — — : 

1 WL iCHT ( I ) = 1 . 0 

NC=-0 

.—C— START-M;^U; GUlOANCi.. FOR FIRST BURN- , ... . 

3 . . IF (NLAST. LO. 1 ) LA 01-1 

IF t DAOSI ( o ) -T1 ) .LL. 2.^0 TYPL { X . >NLaST-1 

T F- C-NLA ST . L u . 1 ) . - 1 O U T = 1- ■ - — — 

lF(NLASr.LU. 1 ) l>T-DAaS( T 1 MLS (t»)~TiMLS(b))/2- 
IFILAST.LO. 1 ) uT=UAiiS I T iMCS( u ) -TI MLS (S ) ) ■ 

CaCC.,.,FORWRDC 0 t ~LT .0 J - - - 

: T iMtiS ( 6 ) 1 -MLS ( o ) -DT , • 

. _ ..... lF(LASr*L.O»l) GG T L' o 


8CB00320 
BCB00330 
—UCU0034 0' 
BCU003BO 
BCl>0u3o 0 
.-LaCiiU-0A70 
BCBOOBbO 
L>CB00390 
...bCBOQAQc. 

BCBU04 1 0 
., bCBOU420 
_Ga.CiiaQ4B.O- 


BCaOOAAO 

oCBOOASO 

tiCaOOAoO 

,..bCu00470 
bCB004<sU 
B.CB 11.049 0 


acBoosoo 

bCBOOolO 

_acuooB2o 


., BCB00S3U 
. L1CB00S40 
-ii CB 0. Q S S O 


O' o 
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FiLii : bCL-i-iJ 




l> 1 


CAMi.iU 1 UC.l: MON I T O /•’ S YST EM 


—M0iVii=0 — - • - - ■ 

CALL CUI OE I U . 0 ) 

IPK1KT=0 

^ CK=-*1*0 

CALL-.CKS^I < lK } 

.C-SL.T-.- INDEX -EOr< -DU TPUi • 

^_I F-C-l-QOT . L<J •- i >— NPU I N I -NPiH f4-T-+-i 

IF ( iUUT.EU . 1 / lAlL NAVoUT < I , NPO i NT} 

C apply COKPECriONb TL* OD* KLiIP '^il (-TiN*UE 00 > AT 
-C-OT— IS -THL ACTUAL H.k>TAT^ WiiicN. WuUlD -.Hz- USi-D 

DO' A I = i ,U 

U0( i ) =u0 ( 1 ) +L . K. •!' ij U { 1 } 

Ui UTtl-)=UU-) - 


BCB00^7. 

BCBOub'3'^ 

^CLi.OODOC 

riCDuOoOw 

DCB0U6 1 L 

ti-Cmj-DB.2X 


UNIT MAGNIT 
FOP. STLLfilNG* 


fOE 


dCdOOgBO 

aCc30U6Lk 

BcaaoBLL 

._...aCB006DC 

tiCb0067o 

_>UVi.OL^ V. 




uMAG=usurt;r (Ouc i ) 

DO 41 1=1*0 

-4^1- U0( 1 )=G0( I )/UMAv. - - - 

_C~APPLY- CUktl^LCT I UNS fu ALL tiMLS LXCLPF THL LAST* SlNcS IT IS 
. C -THE. CLOCK time IN lACK’/A.-^L>S f-iUoL - ... . . . , , 

' ' ■'■■ ■ ■'DO— S— — -•• — — — - — — — . 

S T 1ML£,( I >=T i MLS ( i ) +CK XOT 1MLS( 1 ) 

AM 1 =TK UEMS + VLH ( i » L ) '■!=( T 1 ML S ( o ) ~T I MeS I b ) +T 1 Mir. S ( 4 ) “T J MES ( 


REALLY 


aCSUOSvv 
&CB0070J 
— .jbCBOQ-7.1 ; 
_... li CD 0 0 7B <■> 
.. _.LiCbOG7Dv.' 
43 01.0.0 74 L> 


IT 1MES( I ) ) 

AMO = ( AM0 + A.''5 11^.0 

IF UNLAST *NL . 1 ) uLi Ty D 

— <4C-NC4-1 - 

IFINC.LQ.I) UU lu D.5 


bCBuUTSi.. 
4>) + TIMES ( 2 ) -t>Cb0 0 760 

DCBG0.77C- 

DCiiOOTGG 

.. ....... 1JUJ007S/0 

ucaciuao... 

DCbOOil lu 


..U 


c 


DM I SS=^05QR T ( ( X ( 1 ) -XF ( I ) X ( 2) -Xf (2) )=:'*£+(X(3)~XF(3 

IF( DM I .SS*GT • 1 . . AJiO •Nt . *LL .20 ) OU - TG .Db - 

. lALL GUI uL I u . c ) 

DO t I 1= 1 . o 

GO TO b 

End Or FIRST JURw, MlSSlOiT Ai.iuUNO. 

CALL GUI DL I 0 .0 > .• - • • • - • ■ 

...MODL=j0 

SLT-.FLAG to SIuNAL i.AV 1G/\T UjNAL updatl. 

C RESTUHE T.IL WLIGriTS. 

DO 7 1 =1 *t> 

-7 WEIGHT t i J =GTD( 1) - ^ 

C-CMANGE (JUIljUUNO 4-L>l.‘:;R '••.iSSl^rv D.ACK Tij NORMAL .MODE* 

Di =k<ADS ( •/.!- j I I S) ) 

^ ( T i .'.1 1. s 1 ) “ r i ( 4> ). ) — • — 

rl3-DAaSi T r HES < )-7 ] '4i::.s< 1 > ) 

C2=0AHS( TlMr. S ( b ) -T UU Si 4 ) ) 

E3=i}A(..S ( T I MLSU.'- ) -1 1 2 > > • -- 

. . TT=riMLSt I )-T0 

TO-bi + .fC INT 

TIMES { 2 ) ='J . 0 
TIMES( S)=TO+CL 

TIMES(4)-TTMHS<b )+oL - .. - 

Ti.MESF b) =T IMcS I 4 )+CS 

... . r IMES( G) = 'f I MEb L 4 J +L3 

tt=-ti'mes(6)+tt 




bCUOOS<;u 
... BCBOO.SbC 
bCDOOSAO 
. DCBOOEbu 
.Jb.CdiUu.3 o o 


riCbOOd7C 

aCQ0C33L 

tiCUOOo'^u 

DCbOOyOu 
bLO009lu 
_ ■■_■ c^. CB D. 0-9 E.U 


0CB0 094iC 
bCS0094U 

bCbOOUbO 

tlLSOOyoC 

bCboo97o 

L».Cu 0 0 '5» E u 

bCb00990 

BCUOlOGo 

bCbO.lUl u 

uCuo iuau 
DCbO 1 03u 

._iiCaD,X.0-4.C 

BCD 01 06 t 
6CbO lOoC 

iiCUu.lu7>. 

.. . . bCbUiOE^ 

,6 CE 010 9 C 
acL.a i.T Ju 


BCBCB - 7 


Fim; BCbCB kurtkau 

v-LHC-tw: 1.*:^ } 

AMO = TOTUMS-LU*Vi_H ( 1 , ) 

DU B I =1 

aau )-DTti > 

. .. QO I 1+3 )--Q r t l-i- J ) 

- - XTLMPl =XTl 1 ) 

x-Ti£-MPi: =a.T,< 1-+0 - 

XT( I)-XOt i ) 

. XTC 1+3 ) = -X0 ( H 3 ) 

AiH l ) I - - , 

. 3 XB( I + 3) = -XT£:mp2 

- C-StT BND CUi'.U I T I CiNS • 

^ALL--UV-ALi>-(-X.T-*UU*P i V-*TV-»-l->- 

DU 9 1=1*6 
9 CC ( 1 )=DU( 1 ) 

CALU-C:;iCli - - - 



ucboiistJ 

dCBO 1 1 36 

aciiOi-iAa. 

dCUOl 160 
3CBU 1 1 bO 

tiiiaui. 1 . 7 0 $ 

bCBO 1160 
aCB0il9U 

, - ficuai2 00 

BCUOiAio 

acti0l22u 

LiQLiii.L^L. 

tJCB01240 
bCBO 1260 
BCbOiabO 


. KLTUr^N 

C 

C 3-BURN INHOUNU Miii6lUN. (i-Oi^WAkU MODU), 
C 

C-.SBT UO. Tu UNIT MAyiNl TUOi: * 


UCU01270 

6CB0 1260 

b.Cb U.12J-J 0 

UCBO 1300 
b<-bO 13 I 0 

_UCU0X320 

UCDO 1330 


100 . OMAGsUSOT'^T (Uu < I ) ■;s>’r^+uu I 2 > ❖>{‘2+00’ C 3 ) 

UU-- i-G.i-i-i — : 

101 Q0( 1 )=Q0( 1 )/OMAG 

102 IF( NLAST.liO* 1) LA3T=i 

lF-.<-T IMLISl 2 )-TIMLi6( 1 ) . 1_ lZ . 2 . AD TY36 ( 1 1 )... NLAiT=l .. 

if(nla6t.L'U 1 ) ( 1 i >.K-:6 (2) - ri^flL;s( 1 ) )/2 - 

.. .. lF(LASt .LU. 1 I Ur -M 2 >-ri M6S t 1 J 

CAL-L-i-'URWKlU 0 *U 1 «.l ) . 

I F“C l a 3"t ■^IToTi' ) ■ “ iTJi^TTA“f*^l 

TlMfZS< L jariMBS ( 1 )+LlT 

IF-(LA:>T,i:a*l ) GU TU iOc. . .... 

, N0P=1 

. MUUL = 0 

^CAXU- GUI uci( G..0.) .... 

1 PRINTED 


UCUO 1340 

UCUU1-S5U 

BCQ01360 
bCB01370 
— aoiaiiso 
oCUO 1390 
bCB01400 
^ tu-rioiaif. 

■""B CB 614 2*,, 
6CB0143(, 

BCu.0144l> 

BCBo 1 46u 
B CD 0 I 4 6 U 

aCa0L4 7,0 

bC30 1460 


CK=-1*0 BCB01490 

-CAUL. CSsSLT (UK ) . .......... . BCbOlGOO 

C-ADU UUURL.UTiUivS TU wu, flUuS. bCbOlblO 

- . OU 103 1 = 1.-:. 3CB01620 

uo.(.Xl=iUO ( 1 ) + CAAU..I- ( i - - — - - ^ bOi01i:i3u 

10 3 TlMbSt I ) =T IMLS ( 1 )+;.^T iML3U ) ACK bCoOlS + v 

IP ( KLAST.bU . 1 ,Urs; . UjU I * LU. 1 ) NPU i Nr=NPO 1 NT+ 1 tiCBOlbSu 

IF ( NLA3T .to . 1 • U.i • I'J J 7 ...U • I ) CALL i-.AVUUTt 1 * NPO 1 NT 1 . . . . .. BCUOlbbO 

GU Tu 1 00 , bCaO 1 67o 

106 CALL CbCb BCB01630 

rlt.uhu^ -. ... b.CciO.,i6yo 

£NU bCbOloOO 



BVAL5 - 1 


Subroutine BVAL5 
A. Purpose 

The new BVAL5 subroutine replaces the BVAL5 and BVAL6 subroutines in 
GUIDE 71/6. It calculates the miss in end conditions and partial derivatives 
of the end conditions for either hard or soft constraint missions with up to 
six end condition constraints and free or fixed terminal time. The subroutine 
can also be called (for example, for initializing desired h and e) with 
NBVAL=-1 to calculate the three components of the angular momentum vector h 
and the three components of the eccentricity vector e, pointing toward peri- 
gee with magnitude of eccentricity. BVAL5 calls the subroutine COAST to 
obtain target state XTF at the end of the mission, TIMES(6). 


B. Input/Qutput Definition 


Input Parameter 

Symbol 

Definition 


XF(I) for 1=1 to 3 

r 

Final vehicle 

position 

for 1=4 to 6 

v 

Final vehicle 

velocity 


QF(I) for 1^1 to 3 
for 1=4 to 6 


u Final control vector 

u Final (du/dt) 


PTV(I) for 1=1 to 12 


TV 


K] 


Partial derivatives of T,. with res- 
T T T 'T T^ 

pect to y=(r ,v ,u ,u ) evaluated 
in BUZZ 


Phasing transversality condition 
p(r'*'u)/|r|^ + (v^u) evaluated in BUZZ 


NBVAL 


Flag parameter indicating whether or 
not miss in end conditions and their 
derivatives are to be calculated 


UK 


y Gravitational constant 


C(I) for 1=1 to 3 
for 1=4 to 6 


Desired orbital angular velocity 
e^ Desired eccentricity vector 
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Input Parameter 

Symbol 

Definition 

Z(I,J) 1=1 to 12 


Partial derivatives of final 

J=1 to JMAXl 


y = (r*'^,v'^,u'^,u'*')^ with respect to 
JMAXl Independent variables 

JMAXl 

JMAX 

Number of independent variables,- 

JLAST 

- 

JMAXl + 1 

MODE 

- 

Flag to denote fixed terminal time 
mission 

TIMES(6) 

^f 

Terminal time 

TT 

T 

Target epoch (time at which xj(T) is 
valid) 

XTF(I) for 1=1 to 6 

x-^(T) 

Target state at time T 

WT(I) for 1=1 to 6 

w . 

Diagonal components of weighting matrix 
ranging from 0.0 to 1.0. (W(I)=1.0 if 

the Ith end condition is a hard constra 
W(I)=0.0 if the Ith end condition is un^ 
constrained.) 

Output Parameter 

Symbol 

Definition 

D(I) for 1=1 to 3 

h 

Orbital angular velocity 

for 1=4 to 6 

e 

Eccentricity vector 

DELTC(I) for 1=1 to 6 

Ac ' 

Miss in end conditions 

XTF(I) 

X-j-(tp) 

Target state at tp 

DC{I) for 1=1 to 6 

DC 

Weighted combination of transversal ity 
conditions and misses in end conditions 

E(I,J) for 1=1,6 

3DC‘ 

Partial derivatives of S with respect 

J=l, JMAXl 

. 3 ^ J 

to independent variables 
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C. Method of Computation 

Components of the orbital constants h and v are calculated using the 
expressions 

h = rxv 

^ (rxv)xvi 
11 .- •» 



The subroutine COAST is called to propagate x-j-(T) from T to t^. If a fixed 
terminal time mission is being flown (indicated by M0DE=3), the parameters 
JMAXl and JLAST are each decremented by 1. This has the effect of eliminating - 
the dependent variable corresponding to the change in the transversal ity var- 
iable across the last burn arc. It also has the effect of eliminating terminal 
time as an independent variable and of eliminating the appropriate row and 
column of the E matrix. 

The end condition miss vector ac is composed of scaled components of 
Ah, Ae and Ar lying along the R and K = jfjj vectors and a scaled miss in or- 
bital energy E. 


Here, 



= target ' 

" ^target ‘ ® 


( 2 ) 


and Ac is evaluated at R = r and H = h. This constraint formulation has ex- 
cellent convergence properties for well posed orbit injection and rendezvous 
missions of all geometries. All components of ac are scaled to. have the same 
units as r. 



BVAL5 - 4 


In order to avoid stability problems during the last leg of a mission, 
the problem is formulated so that a weighted combination of fuel use and miss 
in end conditions is minimized. The cost functional 




mldt + 1/2 AC Wac 


(3) 


is minimized. Here W is a 6^6 diagonal weighting matrix and [mj is the rate 
of fuel consumption during burns. Minimizing this cost functional is equi- 
valent to satisfying the costate equations 


Pf 


3AC 
3X , 


V/Ac 


where = (u'*',-u'^) or equivalently the equations 


(I-w)b''^P^ = wac' 


X=x 


(4) 

(5) 


where B is a nonsingular matrix such that 




I 


X=x 


( 6 ) 


and w is a diagonal weighting matrix with ith diagonal component w- related 
to ith diagonal component W^of W by 


A 

1+W, 


(7) 


Whenever an end condition such as phasing is unconstrained, the corresponding 
diagonal component of w is zero. For hard constraints, w=I, the vector B 
is composed of six scaled transversal ity conditions. The sixth component of 
B*^p^ is |rlT^/|h| where T^ is the phasing transversality condition calculated 
in BUZZ. The components of (l-w)B^p^ given in terms of multiplying coefficients 
C.. defined in the code are 

^ i3 
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The DC vector calculated in BVAL5 corresponds to the miss in satisfying Eq. (5) 

DC = wAc - (I-w)b"^P^ (9) 

Partial derivatives of DC with respect to the independent variables c are 
calculated via the chain rule. 


‘aDC) _ f sDCl fix] + fiDC] f^’ 

.3C j “ [ax J [dz] Uq J [a?. 


( 10 ) 


The G” matrix' ih' BVAL'5'W to { DC/ x) negTecting-derivatiyes of - ■ 

scaling factors. From Eq. (8), it can be seen that the second term in Eq. (10) 
is efficiently evaluated by calculating terms such as and mul- 

tiplying by the appropriate C-j coefficients. 
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FORTi-J Ab I--’ I 


•o-'vs;.T.c 4 .) ) 



r-^QNlTOK SYSTEM 




l543-0*b*H4t:i 
C4 l=r( 1 . 0-WT ( 4 ) )/HM 
-C4-^'— t>4-3^KTV- ' ' 

Uo= ( i . O-'.^T ( o ) 

- a F A C = U « W T Y ( 2 . w > i-i J J ■*■ V T U ) . - 


dS2=-WT ( o ) >;^RTVU TU-L)FAC 

BbS=2. O^WT ( S > /'JR '.i-U :>-•?= (CF#i-<TU-i^Tua/K2 ) 
--B64— — 2 • 0-i^w T ( J ) v-'-i V 'v IJ t-UF — 

C32=US^KTV-?CF 


--Cba^-* L*Ua4jR.VU>- 

-bt>-l-=WT-( b ) *RVJ 
Cbl^li ,U-WT(6) )/ ( ) 

-C62-=( i .0-‘^T (b) )4tUR;.’/nM 

CAUCULAIL HARTiALJ L'l- 1)L 11 H RlSPLCT TU R ANO V. 

DO b J 

X( l-.-l )~B I i *0( I ) -Li i vv XU ( j- ) - •- 

L< 1 * i >-*-ci i :r<au( i ) 


G(2,Z)=B214Rt I )+C2i--!:OJl I )-h(C21+RTUD + C22=:'VTU)*UR2^RM^W( U/CF 

_-0t2 ,.iF3) =b??.*V t i ) +^.2/*U t i- ) + l C2 l-’^H'r0O+C22’^VT0> ■1'V( i.) /CF, 

LC J. 1 l=~LiO i^f^VC i ) -D->4-^Rt i )-LbL^?Ut i )+CJ.l*UD( I J 

GIJ • i+0) V I i ) 04 + 

4 04 1 ^0(1 1 + i.i<4-2 ( O.TrLi4j'‘}‘V (-1 ) + C4 I 4VXUO (I.) + C42 VXU (1) 


G<4 , i+3 )s-B43*R( I )-C4 l=n2XOD( I ) -C4 2*«XU ( I J 
G ( b * i ) =8b 1 ( I ) -0 iLj/ ( I ) — Cb I ^00 ( 1 ) -C b2#U ( I ) 

X( b* 1 +3 ) =033 *V (I ) + j I 1 ) +C33’5‘0 ( 1 ) •• 


G(c* * I ) =0 6l’i‘X< I I ) -L‘->i T V U ) 
GCO t l + 3)=-Co i^PrV t i + 3 ) 
_ua-b.-J=l , JMAXI- 


DRU ( J > -R t 1 -J ^ ^ ^ ^ ^ ^ 

Dftjt) ( J )=R ( 1 ) *3 I I 0 , J ) +K ( 2 ) -vZ ( 1 1 . J ) +R( 3 ) *Z ( 1 2 * J ) 

XJVO ( J ) =V t 1 ) *Z ( 7 . J ) +V t ? ) 0 , J ) FV ( 3 ) ^Z( y . 3) ■ ---- 

DMU ( J ) =0 i 1 ) LZ / * J ) 1- J i 2 ) -i« i. U 1 » J ) -»-3 ( 3 ) -^ /. ( J * 3 > 

calculate PAtsTiALL OF UL 1 I <-l RLSHLCf fU CO'STATL .iMi_S 

.^At^T i AL-.-UF CUST/Ff- . « lTJi-R-^5PLUT -TO 1 isiOt-Pt-WOLN T . VARLALiLLS.. 

(FIRST STEP Ur C^^^^1N RULu ) 


£ < i * J > =— Cl 1 +<>'TU ( 3 ) 

- L(2 *3 ) =CL i J i J ) -i-ul.'C’.f UV U t J 1 

L( 3 * 3 ) =C3 1 ^'L'i^UL‘ ( J ) — DRU ( J ) + 1, 33 *0 V U t 3 ) 

E(4 , J ) =C41^ ( Dl 1 ) --l-Z I I b i J ) +D t 2 ( 1 I *3 )+D( 3 J ”/l 1 i- * J ) ) “C4 

= T-Cb i *L)RU J. ( J } - LOS: ^'--DRU t- J ) -Cb3.=:^D VU ( J ) .. 


LVA3U Uil- 


BVAO 1 120 
BVAOI 130 
JiVAZU-lAO^ 


#DHU( J) 


OVAOl ISO 

UVAOl 160 . 

_B.VAA1 LZ Ij- 


BVAOi lao 
3VA01 190 
,*iVA0120 0- 


bVAOlZlO 

iiVA01220 

jiVAcaasa. 


BVA01240 

BVA012b0 

.iiVA0i260- 


._ B YA01270 
BVAOI 260 


b-VAO 1290 


BVA0130U 
BVAO 1310 
JiVA-aX320- 


.. BVAU1330 
,,BVA01340 
_a.VA.aL3bO. 


BVAO 1360 
BVAO 1370 
.BVA 0X360 


BVA01390 

dVAOlAOO 

JX/AOIAIO. 


BVA01420 
BVA01430 
.iJVA,0.X4 4 0 


dVAO I 46u 
BVAO 1460 
_LiVAU.X4u"-U- 


bVAOlASO 
aVAO 1 490 
. aVA 0 1 b 0 o 
BVAOlblO 
BVAO 1520 
ciVAQl U 3Q 


c 

c 

.X 


E ( 6 t J ) = — c 6l*(V( 1 ( ii/,J)l-V( 2>^Z( li ,J)+V(3)*Z( 12*3) 1 C‘2 *D RU ( 3 ) 

ADD U'l PAkVIaL Ur PC 'v^ 7 H RLSPCCT TO STATE TI.‘.1 eS /ARTIAL 

...OF- STATE v; 1 TM- K;t-SPEC7 Tu' I N DEPl-UDLN T VAKlABLiiS- 

D06i-l*O 
DO uR=l*0 

J.) =L(I*3) + ^T1* R ) -rZ-C tv • J.l - - ' 

■ CALCULATE MISS lO Si.iFT C TR A 1 N T S . 

CROSS R 


BVA01b40 

8VAU1S50 

X.VA01S6O 

BVAOlbTO 


BVAOlbSO 

JiVA.O.l-b.9.0 


C _ _ 

X 

X 


I- delta M AuUiviG H 

2. . DLLTA. Lf::-i<CY. 

3* uELTA f.. ALONG' H CRu 

if. DELIA O ALUrX K 

ibCL La • L AL.Eii'vG.. A 




BVAO 1600 
BVAOiOlO 

XVA.O L. 02 . 0 . 
BVAO 1630 


.BVA01640 

_JXVAOi65.U. 




U_- i-. 


BVAL5 - 9 


Flue.: CASiiVj 


r- I '< Ai'- 1 ' i 


<•> i L>(j L SYSTEM 


^ OuLTA.. K,. ALUiNiiS-.-Ji-XHOSS. 

WHti^E Dl_LTA i'i£-P^?[,Sl-N T S UhSlJiEf> MINUS ACTUAL ANU CONSTRAINTS 
ARE SCALt'O TO HAYL UNI i'S UH LENGTH, 


-UELTC C l > = C C ( I ) C 1 ) + CC E- ) C 2.J -t-CC U ) VXX I i>.) /HM 


JiVAUl GO X 
8VA01 67U 
aVA0l6d0 
XV.AO-LG&a, 


-DLUTCC 2 ) =- CCr- ,u^?UK^r;UKn^ C CC -4 ) =»C( 4 )-- 

j, +C t U } -ri_ C _ > -t Lt 6 > ^tCC o ) - i . 0 1/ CC C L I"* C-C 1 J 

a— .-fcC (2 C L i + C(-U-J ^Cl- 3 ) U /.UK^ 


UVA0X700. 

aVAQUZ.iO 

[aVAlil_Z2.a 


DELTC(S)=(C(4;— 0(4J ) *XN( l >+( CCu) — O (S))*XX(2)+(C{6)— U(6 i)#XK(3) BVAO 1730 

DClTCC 4 >-C C t i 1 ) + u( 2 >4R(2 ) + C C 3 )4 k (3 ) ) /HM 3VA01T40 

OEL-TC (G) = CC(4)'-ul4))=!^NC i)-frCC(5) —Li ( G ) ) ( 2 i + C C { 6 ) -D ( G ) ) 1 BVAO 1.750 . 

,...LiriUTC(0) = lXTKl i )^^XK t i ) + X T r t 2 ) t 2 ) +X T F t U ) -I' XK ( o ) UVA0 176 0 

LALCULATu Xcii'oHrLLi COMO li^iATl ONS UF TR AI'JS V LR SAL 1 TY - BVAOl/70 

LUICU A I-iUNEi-- ANL>-...‘4 i LS.. -iN--Liv»J-- CONUl-T I.UNS ••• ' — "" tiVAQ-L7.00_.. 

bVA0l790 
bVAO 1500 

tiVAQia-lO._ 


UC ( I )=-Cl l^HTU+vV r ( I ) =M>EL fC ( 1 ) 

UC C 2 ) =-C 2 1 *KTUD— C 22 ^ 7 TU+W r {2 )^'«DEL TC C 2 ) 

UCC 3 ) = — C 3 t =fckTUD+C 3 .«_=i'NT U-C 33 ^ V TU+A'T I 3 ) *OllT C ( 3 ) 

0 C C 4 ) = — C 4 1 ❖ HTUU + C 4 2 ^. i 1 U + .v T (4 ) ^‘DLL.f C ( 4 J 

DCC 5 ) = C 51 -J'RTUU+CS^J^R i U + C o 3 ^ V TU 4 W F ( b ) * 0 LLT C ( S ) 

- UC ( o ) — Co 1 •>LT. V + ■»' T C-ii ) '‘'L.-L-r c C ^ ) — ■ — • • 


OVA0 1520. 

_.__._BVA0.1&30. 
BV A 01840 


RLTURN 

END 


BVA01850 
13VA01560 

8VA.0,liiZ0._ 

C_<>0<><><><>00<><->n><X0 0<><2<><><>0s><><><><><><><><> OOXXOOOGVAOltibO 

^ ~ ■ .UVA015.90 

--—SU 5 RDUT.INE’ SOLVE lA*LG.*L/) _ , „ ., „ a V A U ,,l, 9 0.Q- 


REAL’J‘3 A( 12 *25 ) .0 

DO 6 N&ltLO 

._M=N+ 1 - 

.. 13iG-N 
.._UU I I^N.Lb 


BVA01910 
BVA01920 
-B.V A 0.1.93 U... 


■^1— ■■!.. 1 (- DA GS'( -"■•A-C L^»i 


1.5.1..0JT.1:: 


. DO 2 JsM,L7 

U = A{ 181G* J )/A( i:.MG .(') > 

_A( ItiIG*J)=ACN,J) 

ACN. J)=0 

AUBIG.N)*ACiJ*NJ 

U i * Lo - — 


bVAO 194 0 

BVA01950 

^^-■^__il VAUT- 9 bQ. 

BVA0197C 
BVAO 19b0 
BVA01.990 


IFC I .EO.N). GO TO !j 
DO 4 K-M.L7 

-A C-I ♦Kl^AC i *K. )“Ai i ,N)^^A(N*K.) 


... 5 -CONTINUE 
.6 CONTINUE 

hlturn - 


. BVA02000 

.._.DVA02010 

LlV A.0 2-0 20.. 

BVA02030 

OVA02040 

UVA02050. . 

_ BVA02060 

BVA02070 

LVAUZ-OO U_ 


END 


&VA02090 
BVA02100 

C 0 <>*^X ><><><><>0 ^><.>^>-^>*..><>0 <><><-> <> O <> O <><><> <><><><> O 0.0 — _iiVA0.2 l-l 0 

^ BVA02 120 

SUUROOTiUL 3VAu4 a'i- » -i r » P F V * X V ♦ MOV AL 1 BVAOLlSu 

— —-IMPLICIT Ri_AL*o ( A-M — : a.VAOZX.40 

C THIS IS A FOUk^CuN -» TR A I N r VLRSlLiN UF fWALG* MISSION BVA02X5O 


IS TO ACHIEVE AN UUU I F WITH GIVEN VALUliS OF SEMlMAJOk AXIS, 
-ECCtNTRlC 1 TY -INCLl .^A 5 1 UN, -AND -AKCU,.^LNT UF. PLRIGEL.- THL.QRBI-TAL 


BVA 02160 

BVAO.2.17.0. 


const ANTS .WHi CH ARE 1 i i A NS.M 1 I Fi-U IN C. IN THE CUMMUN HLOCrC GlOlN- „.,tJVA021SO 

C— .-ARE... MAGN I TUUE- AND CUMPONE: N f OF- ORUiTAL A NG ULAR VLLO v 1 T Y H. -,B VA021 9 0 

-THL— T-H-IRD CGMPGN'EnF i.Ji' A -VECTUP L-PUlN-flNG TO W AP O • P ER I CEN- i.R.— W LTE4— — U VA022D-Ci- 


fill:: CASBvJ hiU<ir<A.; --‘X 



8VAL5 - 10 


CAMbRlDoL MONITUA: SYSTEM 


/CPi }YS/ n-j HU 0, A HO(i * OM t G A * U3LA Tt' 


JiMAuZiLU 

Sy^AO?22< 




T ULL. 

COMMOh 

COMMON /ClNDfiX/ NAf-iu * I APC, JMaX » JM . JMAX 1 * JLAST ,NO •Nl)P*N-- KGOS QVA02:^J^ 

-COMMON /GIO iN/ -XT (<> ) *.r t ,;<0 LU) , To» Ai.lO * VLH.( 1 O * 7. ) * Q 0 ( to ) • CX£U-tlVA0a2A 0 

COMMON /GlUUU r/ LjUi. L o ) • UT lMt:s(6 > ,£( i 2 * 13 ) .ZZ( 1 2 ) ,Z ( 12 • 12 ) . D{ 6) • OVA022toO 

I DUMM(A).SM ' _ -ti\/A02260 

— U-l34L-N:ilO»4 -Xr {-to.) *Ni- ( tol *P-T-\4-{-l2 ) *G(-to.) li.VAUZLZ7U 

1 > CALCULATE D(i) FOR i = I TO 4 4; * A t>VA 02 23 O 


R2=l ,0/( Xr ( 1 ) ;«XF ( 1 )+XF { 2 )^XF (2 ) +,XF ( 3 )«XP (3 ) ) 

)**0-to . , 

RTV = XF ( 1 ) *\r ( A ) + XF I ^' ) ^''XK ( to ) + XF ( J ) '■i'XI- (6 ) 

V2=aF ( A )^rXF ( 4 ) + Xi' I j } ;=Xi - i 3 ) +XF( 6 ) t to) 

AO- Li . ) ^ V 2 / R 2 “U T V .'X X- T V . ) . A-to . to 


C**«< 


O { 2 ) -=^XF ( 1 ) ViXF ( to ) -XF ( > -15 XF ( 4 ) 

H2=D( I )<=U( 1 ) 

.-RTVH-=RTY*i;M 

FV=H2/UK-1 .U/i<M 

0(4 )=NTYK-'>:Xi- { 3 ) hFV !=XF (to ) 

-Ol-4-)^.TV*XF.(3) .A - A I •:. ( o ) / +H-2 # XF (-to) /UK 

V2UR=:V2/UX“RM 

R rvU=H TV/UK 

D-(3)=L<F(3)*V2!Jk~XF ( to J TVU..— — 

IF (NJVAL.Lo.-l ) Rto turn . 

2 > CALCULA1L PARTIAL uLii I V A T I V LS L(I.J) 
>--CALCOLArL- - {Oto4/0N.J •-...{ OC4/toV )•«••■ 


R3R=Xr (3)«ftM-TN2 
V3U=XF (to ) /UK 

■V2H=V2/0{1) 

RTVH-RTV/LX 1 ) 
R2H=1 .0/ CR2'-^l)( i ) ) 
-TRJ^2 • :> ) /.UN ... 


0VA0249O 

bVAO2to0c 

JiVA.U^S.LL 


CS1 = XF t6)*(V2UK+V;?/'JK ) -R TV'S^RtoR 
CS2»2# 0*V3U/R2 

; — F=XF(3 )4<-RM-2.04*N rVU XXF ( to) 

- DU 4 I = 1 • 3 

- G( I )=(=s>XF I 1 +3 ) -H^to l->.At- ( i ) 

4 G(-1-+3)=C52«XF{ I k-toJ+F-A-XFl l ) 

G( 3)-G ( to) +F rVi( 

G( 6 )=G ( 6) +F V 

OU -i J = l.JMAXl ..... 

C***.«C 2.2 > CALCUCA lii ,\*to 

RZ I -XF ( 1 J ( 1 * to > + aF ( : : ) ;-Z ( , j ) +XF ( 3) 3 . j ) 

— RZ-4 “XF (1 )*Z(4>J) + Xr (Z)-v=^(to»U)-l‘AF. (3)-XZ(to*U)..... 

VZ1=XF(4)*Z( 1, J)+XF('to)'^Z(2.3 1+XI-lto)Ac7(3, j) 

VZ4=XF ( 4 ) ^Z{ 4 , J } -(-XF ( ■;. ) ( to , J )-i-XF ( to ) ?=Z( Q * J ) 

-C*-* » ♦ •-< 2.3 > F I i J I ;>H c ALC’JL AT 1 0» 4 G F . l. ( I » J , 

L ( 2 * U J =XF ( to ) ( 1 . J ) -aF ( '( > ->:Z ( 2 , J ) -XF ( 2 ) 7 Z ( 4 . J J +XF ( 1 ) ( 

L( i . 3 ) =V2H-^Uzl -r<T VH VZ i ♦•i{Z4 )+;<2H4 sV/4 

L(3*J) -V2UK:itZ( 3, J )+:<3.<^NZ-i.+.TKto >VZ4-V3U4'< '/Z14RZ4 ) -R T VU*..L( 


bVA02U2 t 
BVA0263 to 

bVA023.4 to 
BVAU2totoO 
OVA02to6G 
-UVA.0Z5.7.U 


BVA02580 
6VA02S90 
avAoatooo 
'JVA02610 
BVA02fa20 
-b-VAO 2Xa30. 


J) 


bVAu264u 
DVA02toS0 
.-•-•aVAU2to60, 
liVAOZaro 
. - bVA026toto 

LtoVA.Q2xji)0, 

3VA0270to 

bVA027l0 


t(4tJ)=G( 1 i»J)^-^J(2)=!'Z(2fJ)-♦'^J(3);sZ(to*J)FG(4)4tZ(4*J^ 

I G < to > =^Z ( to , J ) f to ( o )^Z( to , J ) 

(S»3)-fOF(S)--.S 4( I »J)— Cr (4):^ZC2*J) -toil- (2)’-l‘ZC4*J)+0F( l)=>Z(b * J ) . _-_ii.VA02720 

1 +XF ( 5) =J=Z(7. J ) -XF (4 J (o , J) -XF ( 2 ) *Z ( 1 0 . J ) +XF ( 1 )^Z ( 1 1 , u) . „ JVAO2730 

—SUM-0. 0 -LiVA02740 

— OO— 2 ........ : ^ — LVAu2JLa 0 - 


BVA02290 

3.VAO230U 


C1VA0231 to 
bVA0232to 
^.VA0233to 


BVA02340 

6VA0235U 

-bVA023oto 


bVA02370 
.bVAOZoatot 
JJVA 02-39 0 


BVA02400 

OVA02410 

_aVA02420 


UVA02430 
,'3|t ♦ » V A 0 2 4 4 0 

NtotoUA-024SU 
aVA024toO 
BVA0247O 
-aVAOii.43 0 


bVAL5 - 11 


file: aUki:JVJ 1 CAM'iJt^ iOGii MUNlTQ SYSTEM., 

2 J(lLiA-t>yj^A+i^JV i 'r:, ) << ■ — - — . LVAu^TbO- 

„ 1 ||//(6t J) =SUM BVA0i:77O 

C**^<l/i > CALCULATE MlS:; 1 1\ w-NU CDNDITi(5NS DtL C «**BVA 02? 7b0 

■—IjlpO -A l=l ■■ : aVA027-9^- 

3--ifL*( 1 ,.JS )=C(1 )“U ( I ) ■ ■ - - BVA02S0 0- 

l-i )---T V - •- ^ BVA02iJ10.. 

i ;--XF 1-4 )*CF-( 2 ) -XF.(.2)^tQr.C-4)+XF.( I.) ^yF(B-).)— tLVALlZB^-O- 

#/. return avAoaa^o 

/ i-MH SVAU2B40 





